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Abstract 
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An  ord©r-of -magnitude  analysis  of  the  complete  linefrlzed 
equations  for  a three-dimensional  disturbance  superimposed 
on  a two-dimensicnel  boundary  layer  in  a compressible  Cluid 
shows  that  they  can  be  reduced  to  much  simpler  equations  i.n  the 
first  approximation.,  The  latter  equations  can  then  bo  trans= 
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slonal  disturbance,  so  that  the  mathematical  theory  la  simpli- 
fied considerably,.  In  contrast  to  the  situation  for  an  in- 
comprosalble  fluid,  however,  tl.o  transformation  gives  no  Infor- 
motlon  directly  re^'ardlng  the  relative  importance  of  two-and 
three-dimensional  die turbances , 

First  approximations  to  the  viscous  solutions  of  the 
disturbance  equations  have  been  cbt'tlned  which  are  u:\2f orr.ily 
accui'&te  over  the  whole  boundary  layer*  In  terms  of  these 
new  solutions,  the  relation  for  the  ch&.i*actor Istlc  val\.ies  again 
takes  on  s form  very  much  like  that  in  the  previous  theory  of 
Lees  and  Lin.  For  small  values  of  tlie  non-dimensional  weve- 
apeed  c,  the  pre^iont  results  when  apeclalized  to  two-dimensional 
disturbances  rSviucs  to  those  of  the  earlier  theory,  oo  that 
the  validity  of  th€»  latter  in  such  cases  is  veri-vied. 

Thr«'0-dlnen.-’T onal  disturbances  are  found  to  be  of  little 
Importance  at  low  subsonic  froe-stroam  Mach  numbers,  since  they 
usually  have  a higher  minimum  critical  Reynolds  number  tiian  two- 
dimensional  onec,  except  possibly  under  conditions  of  extreme 
surface  cooling.  As  the  Mach  number  Incr-eases,  three-dlirensicnal 
disturbances  bacoiTio  signi .^icant  under  conditions  which  are  less 
and  less  extreiBb,  until  finally  at  a Mach  number  between  one  and 
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two  they  begin  to  play  the  leading  role  in  many  problems  of 
practical  interests 

At  supersonic  free-stream  Mach  numbers  the  boundary  layer 
can  never  be  completely  stabilized  with  respect  to  rll  three- 
dimensional  disturbances.  There  Is  always  present  a class  of 
unstable  three-dimenalonel  disturbances  no  matter  how  low  the 
temperature  of  the  solid  surface.  However#  for  Mach  nvunbsrs 
up  to  about  two,  aurface  cooling  is  still  a very  effective 
means  of  atabilizlng  the  boundary  layer-  At  higher  Mach  num- 
bers the  amount  of  cooling  uecessery  to  raise  the  minimum  cri- 
tical Reynolds  nianber  significantly  appears  to  become  practiscnlly 
prohibitive. 

Numerical  results  illustrating  these  conclusions  are 
prosanted  for  the  case  of  tha  boundary  layer  on  a flat  platso 
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In  recent  years  there  has  been  consiojrable  interest  in  the 
theory  of  boundary-isyc-r  stability  for  com^  'essiblc  fluids.  Most 
of  the  i*ecsnt  vjorn  on  the  subject  has  been  ased  on  the  original 
tiieoreticol  investigations  of  Leea  and  Tin  (d]  and  [2]).  As  has  oeen 
X-'olntod  cut  on  several  occasions,  there  are  ct.rtein  inadequacies  in 
this  basic  t eory  which  limit  it?  significence  r.nd  wnlch  In  same 
cases  maKe  unreliable  the  results  of  numerical  <-•  Iculstlons  based 
on  it.  -ihe  purpose  of  the  present  investigation  ig  to  clarify  some 
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stability  theory  for  tne  laminar  ooundary  layer  i a compressible 
fluid. 

One  of  the  chief  limitations  of  the  theory  up  o now  has  been 
considers tlon  of  only  two-dimensional  disturbances,  by  analogy 
with  the  procedure  for  an  incompressible  fluid.  Squl,  v [3]  has 
shown  by  means  of  a transform.' tl on  of  the  equations  c the  problem 
that  for  an  incompressible  flviid  every  throe-dlmensiori  l distur- 
bance la  equivalent  to  a two-dimensional  one  at  a lowe;  Reynolds 
numver.  therefore,  for  an  Incompressible  fluid  the  minium  criti- 
cal* Reynolds  number  for  two-dimensional  disturbances  is  he  rainl- 
rnuiTi  critical  Reyrsolds  nuirfoex'  for*  all  dls iur-bancea,  so  tha  tiie 
restriction  to  two-dimensional  disturbances  is  Justified.  If  an  r 
at  empt  is  made  to  prove  a similar  result  for  a c empress lb  e fluid, 
it  Is  found  that  the  equet’ona  obtained  by  ma  king  the  usual  approxi- 
mations (see  [1])  are  too  complicated  to  permit  a tr-ansfor-ma  ion  of 
the  type  used  by  Scuire.  In  the  present  work  an  order-of-magn  tude 
cnalysla  of  the  terms  In  the  complete  equations  for  a thro^-ciinen- 
sicual  disturbance  is  carried  out,  A consistent  approximation 
scheme  based  on  this  now  yields  simplified  equations  of  the  same 
order  of  accuracy  as  the  usual  ones,  but  with  fewer  terms,  so  that 
a traiisformstion  to  equations  of  the  same  form  as  for  a two-dimen- 
sional disturbance  is  possible.  However,  these  traixaformed  two-di- 
me:.slonsl  equations  are  not  the  equetionsof  a proper  two-dimensional 
diaturbaiice,  so  that  no  direct  conclusion  regarding  the  Importance 
of  three-dimensional  disturbances  can  be  made.  For  a compressible 
fluid  the  prciborties  of  three  dimensional  disturbances  have  to  be 
found  by  means  of  a more  elaborate  investigation  than  for  an  Ineom- 
pressiblo  fluid,  'fbo  tr-ansformed  twc-dimenslcnel  equations  lead  to 
a considerable  slmpllf icetlon  mathematically,  however. 
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It  is  li'ound  In  the  present  inveatip;ation  that  for  subsonic 
free-stref,xn  Mach  numbers  thrae-dimenslcnal  disturbances  ere  not 
of  much  Importianco,  but  for  supersonic  Mach  numbers  they  may  pley 
a lecalniy  role  In  certain  clrcurastencos*  In  particular,  for  super- 
sonic Mech  numbers  there  Is  alvays  present  a class  of  unstable 
;hree -dimensional  disturbances  nc  matter  how  low  the  temperature 
cf  the  5clJ.d  surface  is,  sc  that  «t.«h<ii9'ation  of  the 

boundary  layer  In  such  cases,  as  predicted  In  [2],  la  impossible. 

It  Ja  known  that  at  high  Mach  numbers  the  usual  form  of  the 
boundary  layer  theory  becomes  inaccurate  end  [5])»  so  that 

there  Is  some  question  as  to  the  proper  form  of  the  stability 
theory  In  such  cases.  An  investigation  of  this  pi’oblem  [6j  has 
shown  that  as  long  aa  the  proper  velocity  and  temperature  profiles 
for  the  mean  flow  are  used,  the  simplified  stability  equations 
derived  in  this  report  are  still  valid  as  a first  approxlmatisn  at 
high  Mach  numbers.  However,  In  such  cases  the  accuracy  of  a fix*st 
approximation  may  not  be  very  greet,  so  that  higher  approximations 
and  thus  more  complicated  equrrtioiis  may  have  to  be  ooiialdersd  • 

Various  objections  to  sesna  of  the  approximations  made  in  cal- 
cula  tlons  based  on  the  theoretical  developsnerts  of  II]  have  been 
raised.  In  fact,  it  has  been  pointed  out  In  i?J  that  the  original 
theory  of  ll3  has  sometimes  been  used  In  situations  in  which  It 
is  not  valid.  In  particular,  the  viscous  solutions  of  il]  are 
accurate  first  cpproxlma tlons  only  in  the  immediate  neighbourhood 
of  the  point  where  w = c.  However,  at  high  Mach  numbers  o is  often 
quite  large  in  the  atablllt;  calculations,  so  that  in  order  to 
satisfy  conditions  at  the  solid  boundary  tho  viscous  solutions  have 
to  be  accurate  far  away  from  this  point.  For  an  incompressible 
fluid  Tollmlen  IS]  has  obtained  viscous  solutions  which  are 
formly  accurate  all  across  the  boundary  layer.  This  liork  has  been 
extended  to  the  case  of  a compressible  fluid,  and  a new  mathematical 
formulation  of  the  stability  problem  la  terms  of  the  new  viscous 
solutions  has  been  obtained.  In  the  present  report  the  modified 
formulas  for  ths  calculation  of  the  characteristic  values  are  given. 
The  detailed  methematicsl  discusdion  will  be  presented  elsewhere  . 

an  example  of  the  increase  in  complication, in  a second  approxima- 
tion tejrms  must  be  reteinud  whose  coefficients  involve  the  v-component 
of  the  basic  flow  velocity, so  that  to  this  order  of  approximation  the 
flow  cannot  be  considered  to  be  parallel.  This  conclusion  was  reached 
Independently  in  the  thesis  {63  and  in  an  invest Igation  by  Cheng  I 16]. 
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2,  Seiverel  formulation  of  problem 

The  basic  differential  equations  for  the  motion  of  a compressi- 
ble fluid,  in  Cartesirn  tensor  notation,  ere  as  follows: 

Momentum  equations 


t 0%.  \ 


IfL* 

djt*  ' d.^ 
Continuity  equation 
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I = -L 
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i-  i-  ~ y r — u 
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Energy  equation 


(P.3^ 
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ar' 


dJt* 

Equation  of  state 


^ = 9 [i*u 

3«»  / 


+ % e,-^ 


(2.4)  )t>'  = T' 

In  these  equations  e^ ^ is  the  rate-of-s train  tensor  and 
is  the  stress  tensor,  defined  in  the  following  way 


(2.5) 


± f ^ 

2 ^ dJXi 


a-jT// 


(2,6)  T:.  = t cT/y 


All  quantities  ere  dimensional,  as  defined  In  the  lilst  of  Rymbole. 
The  coefficients  of  viscosity  yu*  and  X*,  the  coefficient  of  heat 
conductivity  k*,  the  specific  heats  Cp  and  C^,and  the  Prendtl 
number  xt  are  considered  to  be  functions  of  the  temperature  alone* 
The  gravitational  force  terms  are  omitted,  since  they  can  be  shov^n 
to  be  negligible  for  the  largo  values  of  the  Pr'oude  number  to  be 
expected  in  moat  practical  problems  (see  [l3  )• 

We  consider  a time-independent  basic  flow  end  a small  time  - 
dependent  disturbance.  Thus,  any  quantity  Q(x^,t^)  may  be  wrltter 
as  the  sum  of  a part  and  a part  Q*  (x^, , t^' ) j 
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(2.7 ) Q(^.^  4-  '(yZi  , J:  ^) 

We  are  Interested  in  determinin';  whether  the  distui’bancea  Q’  are 
stable  or  unstable,  that  is,  whether  they  daorerae  or  increaeie  with 
time. Both  the  quantities  Q in  the  disturbed  motion  and  ^ In  the 
steady  motion  satisfy  equations  (2«l)to  (2.4)»  Thus,  if  expressions 
of  the  i'orm  (2,7)  for  all  the  variables  ere  inserted  into  these 
equations  and  quadratic  terms  in  the  small  disturbance  quentltios 
are  neglected,  the  following  system  of  linearized  equrtions  for  the 
fluctuations  Q''  is  obtained: 

(2.8)  li.  5:5  f>^\  ^ MJ.  ^ '^v 


dye. 


(2.9)  In:  ^ I,  (Jli  jl:  pA  ^ 

djt*  d.Z.L\  * J ~ 


(2.10) 
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The  quant5.ties  e 
(2.12) 
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are  given  by 
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Jiiquatlona  {2,(1)  to  form  o systcra  of  lineer  pertlel  differ- 

ential equations  of  the  second  order  for  the  fluctuat"!  ons.  The 
coeff  * cieiits  are  knowri  fuiiotioiiS  glvoti  by  the  steady  flovi. 

3.  Stability  of  parallel  and  nearly  parallel  flows 

When  the  fluid  is  inccxnpressible,  the  basic  steady  flow  cun 
hfl  fitly  parallel;  as  in  the  ceae  of  two  dimensional  flow  through 
a chann»il.  The  formulation  of  the  stability  problei.i  is  then  very 
siruple.  For  compressible  fluius,  however,  these  is  no  such  simple 
basic  flow.  Flows  satisfying  the  boundary-layer  epproxlmf tlon  are 
the  closest  to  parallel  flows.  It  will  be  shown  that  for  such  flows, 
th©  stability  problem  is  essentially  the  sane  as  tha for  a parallel 
flow  which  agrees  locally  with  the  bas’c  flow  of  the  boundary- 
layer  type. 

For  proving  the  theorem  in  Section  4»  ®nd  for  application  to 
boundary  layers  in  three-dlmendao*!  flow,  the  development  will  be 
made  in  a somewhat  more  general  manner  than  usual,  Ias  consider  a 

3t9irii.*  Xlix  Xii  I 1X0  1/  Xi« 

stream  has  components  of  velocity  and  parallel  to  the  plate, 
end  the  flow  Inside  the  boundary  layer  is  parallel  to  the  fi’ee 
stream  at  erch  location  (Xj^,x^).  Indeed,  the  boundary  layer  can  be 
described  in  terms  of  strips  parallel  to  the  direction  of  the 
resultant  free  stream.  Each  strip  behaves  like  a flat  plate  with 
the  loading  edge  normal  to  the  resulting  stream.  The  location  of  the 
leading  edge  coincides  with  that  of  the  particular  strip  in  question. 
Thus,  if  8 is  the  diatunco  of  a point  from  the  leading  edge  in  the 
direction  of  the  free-streara,  the  boundary-layer  thickness  S at 
that  point  is  given  by 

R' 


(3ol) 

where 
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(3.3) 


*•  >&“ 


and  la  the  kiiiematic  vlsicoslty  In  the  free  stream.  As  Is  wall- 

known,  the  quantities  ^ associated  with  the  basic  flow  satisfy  the 
following  order-of-i»iqti^itude  relations  (if  the  free-stream  Mrch 


number 

is  not 

too  large): 

(3.4) 

.(^j  0.  z o(  r;^') 

--OCR) 

0 a.  j, 

\ o'K,  d'x,  / 

(3.5?a) 

{ 

01-  Mj ) ; Mi  - 0 ( Ri^ , 

1 ^ OCR) 

(3.'?6) 

= 0(T.^}  , 

, r = OCF.*) 

Ve  shall  now  simplify  the  equations  for  the  small  fluctuations 
by  means  of  analogous  considerations  regarding  the  fluctuations 
end  their  derivatives.  From  j>hysical  considerations,  and  also  from 
iiifilrhtsiTis cl*  bcu3nd£2*y“vsliis  ssscclr 'ted 

with  the  system  of  disturbance  equations  (i^.6)  to  (2.11),  It  may 
bs  expected  that  the  fluetuations  have  the  typical  boundary-layer 
type  behaviour  in  a region  adjacent  to  the  solid  boundary.  That 
1s,  near  to  the  boundary  viscous  effects  and  therefore  the  highest 
order  derivati’/ea  in  the  direction  normal  to  the  solid  surface  ars 
important,  and  changes  in  the  normal  direction  are  much  more  rapid 
than  in  any  direction  parallel  to  the  surface. 


general  length  scale  of  the  order  of  the  boundary-layer  thickness, 
so  that 


(3.6) 


J.  ^ ^ o(r) 

Q'  d^,  Q d/3t, 


In  the  viscous  region  next  to  the  aurfree,  the  characteristic 
length  in  the  direction  (normal  direction)  is  much  smaller,  and 
is  obtained  by  considerations  anElo^ous  to  the  usual  boundary-layei^ 
arguments  for  the  basic  flow.  It  is  clear  from  (£.6)  that  the 
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equation  of  motion  for  ui  conttlns  a term  of  the  form  ' 

“ f i ^ 

p'  Uj  du^/axj  from  the  Inertial  terms  and  a term  of  thei  fowi 
^*3*uj^/8x2  from  the  viscous  terms.  Thus*  If  these  torms  are  to 
be  of  the  sane  order,  X2-derlvativoa  must  be  lorajer  than  x^-derlva- 
tives.  Now, 

'.3.7.  g = 


•0.  r gi  = 


/ .T*- 


where  R*  la  a factor  indicating  the  difference  of  order  of  magni- 
tude between  and  x^  derivativos! 


(3*9) 


:/( 
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A comparison  of  (3«7)  and  (3.8)  gives  n = 1/2.  The  relation  (3«9) 
Is  now  comparable  to  (3»4)  with  R„  replaced  by  R . . The  relation 

S 0 

analogous  to  (3*5®)  is 

(3.10)  ( u',  o.  ucj  ) • Vi  = 0( 


An  examination  of  the  equations  also  gives  the  following  relations 
(which  are  used  in  Section  4)^ 


(3.11) 


= 0 


^ u-  ) fit  \ u, ) 


Therefore  wo  have  a boundary-layer  phenomenon  associated 

with  the  fluctuations  within  the  boundary  le:/er  for  the  mean  flow. 

In  a thin  layer  next  to  the  solid  surface  with  thickness  of  the 
••1/2 

order  of  d R*’  ^ the  fluctuations  hrve  the  viscous  behavJoi?r 
characterized  by  rapid  variation  in  the  Xg-dlrectlon,  while  out- 
side this  layer  they  ere  usually  Independent  of  viscosity  since 
the  viscous  terms  in  the  equations  are  usually  unimportant  there. 
In  general,  the  order-of -magnitude  relations  (3»9)  «nd  (3»10)  sre 
valid  only  in  the  viscous  layer  at  the  solid  surface  , while  out- 
side it  the  velocity  fluctuations  are  all  of  the  same  order  of 
raacmitude  and  the  characteristic  length  in  all  directions  la  the 


same  (that  la,  of  the  order  of  the  boundary  layer  thickness  (i). 

1/? 

In  all  of  these  order-of-ma  ^nltude  03timatv:.a  C(R  '~)  may  denote 

1/2 

a quantity  not  actually  of  the  order  of  r.ar^nltudG  of  R ' In  the 
strict  sense,  but  smaller.  This  happens  when  certain  quantities 
involved  in  the  comparison  become  zero.  Such  a situation  occurs 
in  certain  layers  of  the  fluid,  where  It  Is  shown  by  a more  com- 
plete analysis  (see  [6])  that  Is  replaced  by  0(R^'^^).  In 

fact,  when  the  disturbance  is  analysed  Into  periodic  components 
(see  Section  5)  the  ordinary  differential  equations  for  the 
amplitudes  of  those  components  examined.  It  is  found  that  such  a 
layer  ("critical  layer")  always  occui'S  at  the  x^-location  in  the 
boundary  layer  where  the  wa/e-apeed  In  the  x^-direction  of  the 
particular-  component  considered  is  equal  to  the  component  of  the 
velocity  of  the  basic  flov-  In  that  directs cn.  There  is  the  possi- 
bility that  in  certain  cases  the  critical  layer  (where  viscosity  is 
also  always  important)  may  occur  outside  the  viscous  layer  at  the 
solid  surface  described  above.  la  such  cases,  relation  (3»9)  with 
0(R^/^)  replaced  by  will  also  be  valid  in  the  critical 

leyer,  so  that  variai'.lons  through  this  layer  In  the  x^-directlon 
will  also  be  very  rspid,  just  os  In  the  layer  at  the  wail.  Uowevor, 
the  order-of -magnitude  analysis  of  the  disturbance  equations  carries 
through  In  just  the  same  way,  sc  that  It  Is  not  necessary  for  pre- 
sent purposes  to  consider  the  more  complicated  properties  of  the 
solutions  of  these  equations. 

l\..  Simplified  equations  for  small  disturbances  in  flows  of  the 

By  making  use  of  the  boundary-layer  approximotloris  outlined 
above  (corresponding  to  relations  ( 3-4) . { 3.5)» ( 3 *9 ) » ( 3*10  ) , ( 3. 11 ) ) 

It  can  be  shown  that  the  equations  for  the  small  disturbances  can 
be  simplified  to  the  following: 
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where 


(4.7)  4i  r 
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Equation  (4.5)  ^-®n  also  be  written  in  the  form 

(I, .6)  r(  !!’■  . LT*\  ir,  ?I"  ) 

' 9jf*  3/X,  d/Z,.  d/Zi  / 

^P\^t*  d4t,  d^sJ  ^ 

where  «r  = C..  Ja‘ /)i  ',  "^hese  equations  are  correct  except  for  terms 

A few  typical  steps  in  the  process  of  simplif lealitjn  are  dis- 
cus 3 ed  below: 

Example  1.  In  the  inertial  terms  of  the  equations  of  motion 
we  have  torms  of  the  kind 

M,  ?5'  -H  Ji^ 

9^.  TP^r. 


'C.% 


-Ha- 


Only  the  first  term  is  Ificluded  in  This  is  justified  by 

the  fact  that 


-1  /2 
which  is  0(R  ) outside  the  viscous  layer  and  C(R  ) inside  the 

viscous  layer,  accordin'?  to  (3»U)»  (3*5)»  (3»9),  (3*10). 

Kxample  £.  None  of  the  terms  arising?  from  the  dlssl^jetlon 

function  ( *^ij  ere  Included  tn  (4.'^).  This  is 

Justified  as  follows.  The  typical  terms  of  the  largest  order  which 

could  be  contributed  to  (4*5)  ®re 


^ o<ti  'd/iCx, 

Their  orders  of  magnitude  can  be  estimated  by  comparing  them  with 
/U*(  is  the  largest  term  in  the  dissipation 

function  for  the  mean  flow.  Thus  > 


i 


{ 


(4.10) 


(4.11) 


2 


On  the  other  hand. 


the  term  p*i^C^dT'‘* 


/ oXj^  retained  on  the  left-hand  side 


( I 


( 


of  be  03  time  ted  as  follows: 


(4.12)  p\  0^  (^|X_  / 

Now,  from  the  energy  equation  for  the  basic  flow,  ”p*S,C  ia  of 

the  same  order  as  the  dissipation  tem^(3U]^/3x2  )®i  Consequently,  the 
various  quantities  given  by  (4*10 ) , ( .'j^.ll ) and  (4«1<^)  are  in  the  ratios 
1:1  :R  outside  the  viscous  layoi*  snu  in  Lne  ratios  ~:h  within 

the  viscous  layer.  We  are  therefore  justified  in  neglecting  the  first 
two  quantities,  which  «ro  the  largest  contributicn  of  tlio  dissipation 


dT 


if.1 


/ At 


function  to  (4- 5). 

In  addition  to  the  approximations  just  described^  we  con  also 
neglect  the  xi  and  dependence  of  the  quantities  in  the  basic  flow, 
which  determine  the  coefficients  in  equations  (4»1)  (i|.6).  This 

spproxlmstion,  which  has  been  repeatedly  discussed  by  vorlous  writers 
on  the  stability  problem,  depends  on  the  boundary-layer  nature  of  the 
meon  flow  and  can  bs  shown  to  be  equivalent  to  nogiaetlng  terms  of 


the  order  of  R"’‘*‘(see  [6]). 

The  system  of  equations  (4*1)  to  (It. 6)  is  essentially  of  the  form 
expected  from  a strictly  parsllel  flow.  One  Important  characteristic  of 
the  system  ia  the  absence  of  dissipative  terms  in  the  energy  equation  . 
This  makes  it  possible,  as  we  shall  see  in  Section  5»to  reduce  the 
problem  for  a three-dimensional  disturbance  to  a two  dimensional  problem 
In  previous  invest!, getlons  such  a reduction  was  not  possible  because 
of  certain  terms  retained  In  the  equations  for  the  disturbances  which 
prevent  a proper  tranaf oz^aticn.  Thes  e terms  arc  now  shown  to  be 
negligible. 

The  discussion  up  to  now  has  been  restricted  to  the  case  of 
moderate  Mach  nunibers,  for  which  relations  (3.1),  (3. 4)  and  (3-?) 
valid,  t’or  very  high  supersonic  Mcch  numbers  the  order-of-m:  gnitude 
relations  in  the  boundary  layer  corresponding  to  (3,1 ) , ( 3.4 ), and( 3, 5) 
are  different  (see  [^]},and  conaeouently  those  for  the  disturbances, 
correspondi:ng  to  (3.9)»(  3«3.C)ar»d  (3*11)  are  different.  An  order-of -mag- 
nitude analysis  (see  [6])  quite  slianar  to  that  described  for  the 
case  of  low  Mach  numbers,  shows  that  the  disturbance  equations  are 
again,  in  the  first  approximation,  of  just  the  same  form  as(4.1)  tc>  (4.6). 

Hie  proper  velocity  and  temperature  distributions  given  by  a modified 
boundar-y-lsyer  theory  (see  (41 ) would  have  to  be  used,  of  course,  at 
very  high  Mcch  numbers.  The  error  estimate  for  the  neglected  terras  in 
the  disturbance  equations  ia  now  different,  however,  being  ex- 
pressed not  In  terms  of  R but  In  ter*ms  of  a Reynolds  number 


-10- 


(4.13)  R 

based  on  fluid  properties  at  some  ropresentctlve  torioerature 
’f,.  Inside  the  bound.' ry  lay^sr.  Since  T^'  la  roughly  proportional 
to  Ht®  at  high  Mach  numbers,  R is  much  less  than  R in  such  cases, 
so  that  the  errors  in  the  simplified  equrtlons  ^which  are 
C( might  not  be  very  small.  A first  approximation  might 
not  be  very  accurate  at  high  Mach  numbers,  therefore . 

The  present  analysis  is  carried  out  for  any  boundary  layer 
over  a flat  surface,  i.hether  there  is  e pressure  gradient  or  not. 
Thus,  the  results  are  valid,  not  only  for  the  zero  pressure- 
gradient  flow  over  a flat  plate,  but  else  for  the  flows  over  wecfje— 
shaped  profiles.  Although  a pressure  gradient  Is  associated  with 
the  latter  flows  at  subsonic  Mach  numbers*  terms  due  to  pressure 
gradient  do  not  appear  explicitly  in  the  disturbance  equations  to 
the  accepted  order  of  approxlme tion. 

Although  the  eimpllfied  equations  (4*i)  to  (4..^)  obtained 
by  rather  crude  order-of-mag;iltude  arguments,  an  a posteriori  check 
verifies  that  the  approximations made  are  valid  in  the  stability 
Inveatlgetlon  of  the  boundary  layer  over  a solid  surface. 

A move  careful  examination  is  neces.-ary  when  there  is  surfaco  curva- 
ture. Preliminary  results  indicate  tbet  for  most  three-dimenslorial 
disturbances  curvature  effects  are  negligible  ; hen  the  curvature 
la  smell,  but  that  for-  those  disturbances  with  propagation  direction 
nearly  noriaal  to  the  frec-stroam  direction  (se„  Section  5)  additional 
terms  have  to  be  rtjtained  in  equeMon  (4»i^). 

Free  boundary  layers  (such  as  wakes  and  Jets)  are  excluded  from  all 
of  the  above  discussion.  The  effect  of  viscosity  on  the  stability  of 
such  flows  Is  quite  different.  Further  investigation  of  these  cases 
Is  being  carried  out. 


! 


I. 
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5.  Periodic  solution. 3 o£  the  equations  fcx-  smell  disturbances 

The  oqunf  ona  (/;.!)  to  (4*7)  contain  3c^,x^  and  t*  only  thr  ju^h 
the  partial  derivatives,  ao  that  solutions  of  the  form 

Q ~ q (fix)  CAP  i (a,  a,  t otj  'Zj  ^ or,  c*J  *)  | 

J 

are  possible.  It  nay  be  expected  that  a distiirbaues  of  a fairly 
gsnersl  nature  ean  ’00  anal^aod  into  normal  modes  of  thJs  type,  with  a-j^ 
and  real  and  c*  in  'general  emplex.  A normal  mode  is  damped  ;i.f 
c*-^  0,  neutr.'ii  if  = u,  ana  sell- excited  If  Cj>  0. 

Evidently  solutions  of  the  form  '5.1)  may  be  regarded  as  waves 
propagating  in  a direction  in  the  (x, ,xo)  plane  specified  by  the 
number* (a^,a^) (with  changing  amplitude  In  time  If  c la  complex). 

Consider  in  particular  the  case  a^=o.  Then,  the  terms  involving 
the  X,- derivatives  drop  out  of  the  equations  (4«l)to  (4»7)  and  (.ho 
oquf.tions  for  Uj,U2>p  #P  »T  are  exactly  the  same  as  if  the  com- 
ponent of  the  basic  flow  were  zero.  Thus,  for  waves  propagatin'^  in 
the  xj-direction  only  the  component  of  the  basic  flow  (besides 
the  distributions  of  the  scalar  quantities  "p*  end  T*)  has  any  iii- 
fluence.  Now,  s rotation  of  the  co-ordinate  system  in  the  (x^jX-j) 
piano  will  not  change  the  system  of  equations  (4.i)to  (4.7).  dence, 
we  nay  conclude  that  for  waves  propagating  in  any  direction  only  the 
componiBint  of  the  basic  flow  in  that  direction  has  to  be  considered.  In 
addition  to  the  distributions  of  the  scalar  quantities.  Every  periodic 
throe -dimensional  disturbance  may  be  treated  In  terms  of  £ correspond- 
ing two-dimensional  pi>oblem.  One  need  only  take  the  pro;>er  comp<>nent 
of  velocity  in  the  basic  flow  to  carry  out  the  simplification. 

If  we  reduce  the  above  equttions  to  dimensionless  forms.  It  la 
clear  that,  the  wovos  propagating  in  the  direction  of  the  resultant 
free-stream  velocity  are  associated  with  the  Reynolds  nuraoer  R end  the 
Mach  number  of  the  basic  flow.  For  a wave  propagating  at  an  engle 
9 relative  to  this  direction,  these  ntjunbers  are  reduced  to 


(5.2) 


R = 


0 

1 « V- 


C?J!  0 


Mf  = M,  cos  9 


A direct  transformation  verifying  th'a  statement  will  be  given  In  the 
next  section. 

The  advantage  of  the  present  tr^-atment  lies  in  its  generrllty  and 
simplicity.  V.hether  the  incomin.g  stream  is  oblique  or  normal  to  the 
leading  edge  of  the  flat  plr  te  makes  no  difference  in  the  discussion. 


Vhen  specialised  to  the  incoinprens Ible  case,  the  present  results 
cover  both  the  original  vor^.  of  3:iulre  l3l  end  the  esse  of  an  oblique 
InccrniTVP;  x'low  discussed  by  Kuethe  [93  • ^or  cvio-dimenslonal  parallel 
flows.  Squire  obtained  a set  of  transformation  formulas  leadln"  to 
a ccnclii  :icn  equivalent  to  Snoplftmentary  diacusalons  are 

required  In  Squire's  method.  Kuethe,  doallnq  with  the  case  of  the 
oblique  incoming?  stream  by  Squire's  method,  noted  the  variation  of 
^effective*'  Reynolds  number  with  different  directions.  He  did  not, 
hoi  ever,  arrive  at  the  conclusion  reached  hero  that  the  maximum 
"effective"  Reynolds  number  is  reached  by  waves  propagating  In  the 
direction  of  the  freo-stream.  Consdqueritlys  these  waves  first  give 
rise  to  instability  as  the  Reynolds  number  of  the  flow  is  increased. 
6.  Equations  for  the  amplitude  functions 

Consider  now  tho  stability  of  a two-dimensional  boundary  layer 
over  a flat  pli to  coinciding  with  the  {x^,x^)  plane,  and  Introduce 
the  dimensionless  variables  given  in  the  List  of  Symbols.  The  free- 
stream  velocity  is  new  U^,  since  = 0,  Tha  free-atreata  reference 
quantities  U^,  ^I'^l'/^l  evaluated  at  a fixed  point 

.Xj^  = x^(that  is,  at  a definite  distance  from  the  leading  edge  of  the 
plate)  In  the  region  of  the  boundary  layer  that  is  of  .interest. 


..1,3- 


Tlio  dimensionless  velocity  and  temperature  distributions,  w - 
and  T = f depend  on  a number  of  factors,  in  particular,  the  boundary 
conditions  at  the  wall  and  the  free-atroam  Mach  number  It,. 

V/e  consider  a small  disturbance  of  the  boundary  layer  in  the  manner 
discussed  previously,  with  solutionis  of  the  form 


{6.1)  Q»(x,y,z,t)  - q(y)  exp  j i(ax  4-  pz  - a ct)j 

for  the  dl,3turbanco  quantities,  as  indicated  in  thje  Llct  of  3ynibcl;;.W4»n 
relations  of  the  typo  (6,1)  are  injeserted  into  equt-tions  (ij.,1)  to  (lj.,7), 
the  following  equations  for  the  amplitude  functions 

f,  a (f,  h,  7f..  r,  8 are  obtained: 


(6.2) 

(6.3) 

(6.4) 

(6-5) 

(6.6) 

(6,7) 
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f,  ia  (w-c)r  + p»a<|>  + p(la  f + ip  h + ) = 0 

pf  ia(w-c)6  -5-  T»a<p|  = -{I'  -l)(ia  f i ph.  a'|*  ) X- 
^ J o-R 


p-'f 


e 


^ a.  ^ 

P P ^ 


The'" combination  of  terms  af  ph  suggests  its  replacement  by  s novi 
variable  a f.  Indeed,  the  following  fcransf ormstions  convert  equations 
(6.2)  to  (6.7)  to  those  of  an  equivalent  two-dimensional  problem: 


(6.3) 


a f + P h = a f 

a <|)  = a ({> 

a T.  - a 'ii 


1 


c = c 

= a 

aR  ~ n R 
aM,  ™ cT  M, 


a o 


The  two-dimensional  equations  are; 


(6«9)  P < itt(w-*o  )f  + VI*  CL  Tr  ^ f” 

( J <fM|  R 

(6.10)  p^ia®(w  - c)^j  » +>^a<j»” 

(Cell)  l5(w  “ f:)r  + p*a9  + p{iaf  + ~ 0 

(6„12)  pjla{v  c)Z  + Tiap)  * -(y-l)(laf  + af?  ) + 

^ ~ f» 

Evidently  the  transTormations  (6.8)  are  equivalent  to  (5-2)  with 
eoa  © = a/\/a*  V p ao  that  © does  in  fact  give  the  direction  of 
propagation  of  the  vraves  (6ol).  The  dimensionless  wave  speed  e ia  not 
changedp  which  Indicates  that  the  phase  velocity  In  the  x-direction  is 
) actually  sec  © times  the  phase  velocity  In  the  direction  of  the  wave 
normal. 


To  complete  the  aiialyalsj  it  is  necessary  to  vGrl.fy  the  30\2l\’'3 ienoe 
cf  the  boundary  conditionso  In  the  stability  problem  the  boundary  con  ■ 
dlticns  are  homo  gene  ouSs  so  that  there  ie  no  difficulty  in  this  r-e,i.poc!;, 
The  homcgeneouo  boundary  cond5.tions  of  the  thrae^diiaensional  problem 
clearly  transform  to  the  proper  homogoneous  boimdary  conditions  cf  the 
squlvalent  two-dimensional  problem.  For  forced  oscillations  however;,  i 
le  obvious  that  the  equivalence  would  nave  to  be  more  carefully  ex&.minod 
We  have  now  shown  that  the  boundary  value  problem  assoc ia ted  with 
the  system  of  eqtiations  (6.2)  to  (6.7)  for  a thi’es^-dimenaional  disturb - 
isnce  Is  equivalent  to  that  anaociated  with  the  system  of  two'-d.imensionni 
0quation.s  (6.9)  to  (6. 13)  Thus  the  stability  problem  for  three- 


dimensional  disturbances  can  always  be  formulated  as  a t'wo-diraensicnal 


problem  which  is  mathematically  the  same  as  that  for  two"dimensional 
disturbances.  It  should  be  realized^,  howeverj  that  the  equivalence  Is 
p\irely  mathematical  since  equations  (6.9)  to  (6.13)  do  not  represent  a 
proper  two-^dlmenslonal  dlsturbanoot  There  la  a decisive  difference  be~ 
tween  the  present  results  and  those  of  Squire  for  anv  incompressible  fluid. 
In  'Che  latter  casop  because  the  nach  number  does  noz  enzevy  ir  is 
possible  to  transform  the  equations  for  a three-o’lmensional  distui’bance 
to  equations  which  actually  do  represent  a real  two“»di.menaional  diatuirb^- 
anee  at  a lower  Reynolds  number*.  It  is  then  possible  to  conclude  that 
three^dimenalcnal  disturbances  can  be  ignored  whan  the  mlriimnw  critical 
Reynolds  number  is  being  investigated.  For  a eonprssslble  fluid  this  is 
not  true  because  of  the  necessity  cf  transforming  the  free»streara  Mach 
number.  It  la  not  even  true  that  every  three-dlmensxonc.1  disturbance  is 


tc  oL  OHS  HiC'iJt  S.  yipd  i ii*  W0l.i 

as  Reynolds  nxjmberf,  because  with  a reduced  Mach  niuaber  the  velocity  and 
temperature  distributions  would  also  have  to  be  modified.  In  the 
equivalence  discussed  above  these  distributions  must  remain  the 

To  apply  the  results  of  a two'-dimensional  analysis  to  the  th’^eo  - 
dimensional  problem  we  may  use  the  transformations  (6.8)0  with  the 

Q asaOciatou  wxlill  ul~i6  tW0''-*d‘iwi6ii6iGT*£.i  prOwiGTiio  Xo  Xi:i  COX*^- 

venlent,  therefore,  to  think  of  three-dimensional  dlaturbiMices  as  waves 
prapagatlnti  in  various  directions.  Per  each  value  of  (that  is  for 
each  value  of  the  direction  angle  of  the  waves)  the  stability  problem 
is  investigated  by  means  of  equations  (6.9)  to  (6d3)  with  the  usual 
methods  of  the  two-dimensional  theory.  Thus^  for  each  value  of  Mj^ 
there  is  a neutral  curve  in  the  (a  R)  plane  which  tranef orms  to  one 


J-O 


in  Lha  (a,R)  plane  by  means  of  the  transformation  (6.0)^  The  family  of 
all  snch  neutral  crirvos  far  all  possible  values  of  in  the  range 
0 < ;<  gives  the  stability  characteristics  of  tbroe^dxmensior.al 

disturbances,. 


7. 


for 


Fortaulatlpn  of  ^he  charaoterlstic^value  problem 
Having  established  a two-dimensional  fonnulation  of  the  problem 
three-dimensional  distvirbanceS;,  wo  may  make  use  if  the  theory  of  Lees 
Lin  [1.1  for  further'  investigations,  fica^everc.  some  modifi cations  are 


necessary  when  the  wave<»8pesd  of  the  disturbance  Is  not  eriiall  ■30B5)ared 
with  the  free-stream  speedy  5'or  supersonic  f roe'^strea®!  volocitloH  this 
la  usually  the  most  importiuit  casOo  Detailed  mathematical  diacusRions  of 
the  neceedary  modifications  will  bs  given  in  a separate  report-  Here  we 
shall  merely  give  a sumznary  of  the  theoretical  argvavents  and  indiijate 
how  the  final  formulas  differ  from  the  existing  ones- 

To  determine  the  analytical  nature  or  the  system  of  equations  (6.9) 


to  (6-13)  it  is  M.n  Introduce  the  varisfc-lss 


(7.1) 
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Then  the  equations  transform  tc  a system  of  .six  equations  of  the  i”irst 
order  - xne  xaever^  wnexi  reuuceu  uo  iicnual  f ona^,  iiave  coeLflcient  o 
analytic  In  the  Independent  variable  y and  the  parameters  M?-  aR  a® 


and  c. 


Six  homogeneous  boundary  conditions  must  be  satisfied  for  a natural 
oscillation.  These  are  given  by 

f 5 ^9  0 bounded  as  y ->•  °°- 

'j  f(<^)  ^ |*(Q)  — C,r , and  eitiioy  |0  jf(!l ) 0 for  an  insulated 

V ’.y 


^ watLl*;  Or  0(0)  “ ‘0  for* 'fixed '.tew^er attire  at  the  , wall. 
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The  condition  of  houndeclnass  for  y to  can  bo  put  into  an 
analytical  form  with  the  aid  of  the  differential  equations  (6,3)  tc 
{6c  13),,  For  large  y viscosity  effects  become  negligible,  so  that  the 
teznns  of  order  l/R  can  be  neglected,,  and  the  equations  reduce  to 
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«>  this  la  easily  seen  to  become 
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Equation  (?»4)  solutions  py)  s where  (for  real  c)  0 
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p ■ ®vt:  5^(1.,,,. 


ThuSs  ^01*  real  c the  oondition  of  boundedness  at  inf  inity  la 
autcsnatioally  satisfied  if 


(7.6) 


c < 1 ~ -i- 
M. 

"1 


In  this  ease  the  characteristic  values  are  continuous^  and  the  solutions 
represent  disturbances  travelling  at  supersonic  speeds  relative  to  the 
soaaponent  of  the  fr©c-=»stream  velocity  in  the  direction  of  propagation 
of  the  disturbances,.  These  supersonic  distxirbancaa  eorreapond  to  sound 
waves  propagating  into  or  out  of  the  boundary  layer  {see  [1]). 

For  subsonic  disturbances  with 
(7.7)  c > 1 -i- 

'tl'w  nirtTiit-.l  Qn  r>f*  ('7.‘X^  wVilftVi  ■?  i»  bmiru^Ad  n't:  Inr-iyil  ■fc-o’  Hm  t*Tnf  ■„  R tt  ^ 

Therefore  the  boundary  cenditlen  for  y ->  >»  16 


(7.8) 
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In  thia  ease' the  characteristic  values  are  discrete,  and  uro  doTjerruinec 


from  8.  secular  ©quEtion  which  can  b©  reduced  to  the  form 


(7.9)  E(S.a.  = F(z) 


In  this  equation. 


(7.10)  pfz)  - 


where 


(7.11)  2 


P(z)  Is  the  function  of  Tletjens 


The  function  E(a,c,M|^*)  is  given  by 
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The  quantities  ?lw»^2w  are  the  wall  values  and  4ll>^21  free~stre 
values  of  two  linearly  independent  solutions  of  equation  (7.3). 

The  principal  differences  between  the  present  formulation  and 
that  given  by  Lees  and  Lin  are  in  the  formula  (7«11)  for  z and  the 
first  factor  in  (7ol2)*  The  results  of  the  previous  theory  are  ob- 
tained by  assuming  c to  be  small  and  keeping  only  the  first  term  in 
each  of  the  power  series  expans  ions  of  w-o  and  V about  y^«  In  thr; 
supersonic  case,  when  c is  not  small,  svich  a procedure  does  not 
represent  a good  approximation  for  our  purposes * 
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^^otermlna tlon  of  the  curves  of  neutral  stability  for  oblique 
waves. 

Consider  oblique  waves  -.ith  ci  recti  on  an."lc  arc  cos 
whoi'e  is  seme  value  in  the  range  0 ^ < Mjj.The  condition  (7.9) 
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Where  z,X,  and  .in  + iv  are  given  by 
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(6,2a) 


(6.2b) 


(8.2c)  u t i y\r 
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4^'  + 6t|ji  <p,. 


For  some  purposes  it  is  more  convenient  to  solve  (8.1)  for 


U + IV, 

so  thar 

(8.3) 
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(8. kb) 

Before  discussing  the  determination  of  the  neutr-el  curve  we 
shall  write  down  the  required  formulas  and  make  certain  trerisformatlons 
which  are  useful  in  the  calculations.  The  quantities  4ll  i>21 
ktX'v  the  f ree— 3 treaiTi 
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volues  of  two  linearly  independent  solutions  of  the.  inviscld  eQuntion 


(G.5)  i.  - ^'1'  1 =0 

ttf]  T-  5 

'Iheae  two  solutions  may  be  obtained  as  power  series  in  a*,  as 
described  in  (l3,  so  that 
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Kor  values  of  c that  are  not  too  largo,  thfe  integrals  (6.?a) 
and  (C.?b)  can  be  put  into  forms  that  aro  more  suitable  for  calcu- 
lation  purposes.  For  exanple,  in  Kp  we  use  the  relation  . 

( r T-M.V-c)*- 


(O.G)  \ 


-n,  i _ \ 

(>v-c)^  1 


to  tr^ris jT ci‘i'1  -^2 
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(8«9)  Kj  = H, K,  Nz  ; vy/here  A4 


Jr  ^ ^ 


Thus  Ng  has  on  integrand 

,eao, 

vhich  is  not  numericslly  large  for  any  value  of  y,  so  that  «p  may  be 
expected  to  be  nuraerically  smaller  than  Kp.  Similar  transfomr ticns 
are  nsoo  ror  tne  other  integrals  as  fdlows; 


= K,  H.,-,  -R,  , -1  , Wj 


(Coll) 
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In  general,  ^ and  M VJill  be  numerically  smaller  than  th«  uorre- 
spondlng  and  for  small  c„  The  preceding  transformations  were 
first  used  in  the  stability  theory  for  incompressible  fluids  [lOj, 
where  c is  usually  fairly  small,  and  then  extended  to  the  or se  of 
compresn ible  fluids  in  [2j.  Obviously  they  are  not  useful  for 
reasonably  large  values  of  c,  in  which  erse  other  trangi'orrn-' tlons 
may  be  moro  appropriate.  This  situation  occurs  for  large  enough  Much 
numbers  - for  greater  than  about  twoj  perhaps.  VJhen  becomes  very 
large.  In  fact,  c is  often  very  close  to  unity’-. 


Relation  (8.2c)  may  now  be  written  as 
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Now  lot  k = ~/*^’/' i~-  suppose  that  kH,a  < l.Then  G(a,c;i.<iy  be 

(I  -6  J. 

expanded  as  a power  series  in  a,  so  that  on  takia^^  I'-eal  end  imaginary 


parts  of 

( o • 12 ) ws 

find 

>a  -- 

y 

«ve  a L 

- A H,  5!  f ( 

(8.15b) 

t 

^ -r  f' 

r ~ 

i r 

rt  . 2. 

"1  'I 

/u  ^ — w 

i-  £ 1 - 2 

t1;.cr  + - 

• • • 

I ' 

a L 

■ i 

[I 

The 

expansion 

of  G 

as  a power 

series  is 

valid 

(8.15) 
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which  is  in  -enerr  i true  if  kH.^  a<  1.  The  latter  inequality  holds 
for  all  cases  at  subsonic  and  slightly  supersonic  Mach  nuiabers,  but 
at  high  supersonic  Mach  numbers  it  appears  to  hold  only  for  oases 
of  extreme  cooling  since  otherwise  k jnsy  be  very  large. 

Tn  the  i>resent  worK  we  shall  neglect  all  terms  of  order 
a®  in  (G.ll+a)  and  (G.l4b),for  the  ss.^e  of  simplicity,  so  that 
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(8. I6e) 
(6.16  b) 
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In  the  cElculatis^na  of  12]  the  same  approxime t Ion  {C.l6b)  for 
V is  used,  but  terms  of  order  a-*  are  retained  in  the  relation 
(S.i.'is)  for  XX.  The  results  show,  however,  that  very  oft«>^  the  error 
corresi^onding  to  the  neglect  of  terms  of  order  in  (8.l6a)  Is 
no  larger  then  t’aal  in  (8.16b).  In  any  everit,  the  approximations 
(8.l6a)  end  (8.l6b)  appear  to  be  reasonably  adequate  for  all  sub- 
sonic and  slightly  supersonic  Mach  numbers,  with  er-rors  within  the 
limits  of  error  of  the  asymptotic  theory.  Moreover,  they  are  very 
accurate  at  any  Mach  number  when  a ^ very  small  (which  occurs,  for 
example,  when  cooling  at  the  solid  surface  Is  sufficient  - see  Section 
10).  For  Mach  numters  greater  thaxi  about  two,  however,  this  is  true 
only  for  cases  of  fjxtreme  cooling.  The  present  approximations,  there- 
fore, are  probably  not  accurate  at  such  high  Mf ch  numbers  in  most 
cases  of  practical  importance,  and  the  method  of  numerical  calculation 
may  have  to  be  revised  . 

Novj  let  us  consider  the  procedure  for  finding  the  oharaeterlsti  c 
values  (a,R,c)  on  the  neutral  curve.  The  various  steps  cen  be 

^ ^ .m-  t-...  — A.’l. 

.d.  ^ .Mm  WW  m»m  ^ ^ ^ J ^ UJ.  Vi/X-iO  f MX  UX  J-  lO 

aid  of  the  condition  for  the  characteristic  values  in  the  form 

\ ^ f A U i-i  \ \J  * X i • 

}*  Cv  oiiC  values  of  c that  occur  In  the  calculations  X(c)  Is  always 
quite  small,- 30  that  i{r^(z,c)and  c ) are  approximated  fairly  closely 

by  ^^(z)end  ]|[^(2).  In  fact,  for  c = 0 , z,0  ) = ^^( z ) and  iji ^( z/) ) = 3],  ( z ) . 


*A11  existing  calculations  suffer  from  limitations  of  this  kind.  In 
some  publlcf tlonss numerical  calculations  were  carried  out  by  using 
methods  adequate  only  for  small  c and  a,  even  though  the  resulting 
values  of  these  quantities  were  not  small.  Such  results  must  be  treste 
with  reserve. 


( I 
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Tlasrefore,  for  each  value  wi’  e,  tha  greph«  of  Ijij-  and  against 
z have  very  much  the  same  shapes  as  those  of  and  respectively 
(see  Figure  1),  The  maximum  ^ijj^(c)  on  the  graph  of  ti(z>o)  against 
z Is  generally  quit©  close  to  the  raexlmura  value  « 0,580  at  z = 3,22 
on  thi  curve. 

For  Maoh  nombers  that  are  not  greater  than  about  tvo,  the  graph 
of  v^(c)  against  c has  one  of  the  typical  forms  shown  In  Figure  2(a) 


/ ^ V\  - _ - - - 

.fcW 


frors  the  data  glveu  ia  i2j)#  In  general^  the  greater 


the  cooling  (that  is,  the  lower  the  wall  temperature)  the  higher  the 
Vq  curve  Is  In  the  (c,Vq)  plane.  For  larger  Mach  numbers  the 
curves  may  have  more  complicated  shapes,  as  indicated  in  Figure  2(b). 
For  the  present,  we  shall  consider  only  Mach  numbers  for  Which  the 
curves  have  the  simple  behaviour  shown  in  Figure  2(a), 

The  pairs  of  values  (a,R)  are  to  be  found  for  a series  of  real 
values  of  ri.  The  mlnliTsum  value  of  c for  subsonic  Mach  mmibara  ip 
c •=  0,  and  for  supersonic  Mach  numbers  is  either  c = 0 or 
c 1—X/K^  ir  > io 

For  each  value  of  c,  relation  (8,3)  is  solved  for  z and  u with 
the  aid  of  (6,l6b).  First,  z is  obtained  by  solving  v (o ) *=  ljF,(£,e ) 
which  Is  ©quivalenito  finding  the  value  of  z corresponding  to  the 
point  of  intersection  of  the  curve  l{r^(z,c)  against  z with  a straight 
line  parallel  to  the  z-axls  at  a distance  v.( c ) from  \ t*  til® 

behaviour  of  the  graph  of  ( z,o  ){whlch  is  similar  to  the  graph 
of  in  Figure  1'  we  see  that  If  v^(c)  < 0 there  is  one 

intersection  point  (one  value  of  z),  and  if  Vg(c)  > 0 there  ere 
ganerally  two  intersection  points  (two  values  of  z)  if  Vq(c)  la 
not  too  large.  For  sorae  value  of  c lerga  enough  it  is  deer  the  t - 
the  straight  line  will  he  lust  tang?^nt  to  the  corresponding 
graph  of  T^i(zp)  at  it?  moxlmum  point,  and  for  larger  values  of  c 
there  will  be  no  intersection  points  at  all.  This  maximum  value 
of  therefore,  i;j  given  by  the  lutsrsectloii  of  the  graph  of 
with  that  of  v^(c)(see  Figure  2).  Gorrespondlng  to  it 
there  is  Just  one  value  of  z,  Vltli  the  value  of  z 


t 


deternii.neu,  u = ^^(z,c)  is  obtained  frori  the  graph  of 
(which  l3  similar  to  the  grs ph  of  in  Figure  1).  rh«n  a 

follows  by  solving  (8.l6a),  since  all  the  quantities  in  this 
equation  besides  a are  now  Icnown  functions  of  c,  and  finally  R 
is  obtained  from  (6.2a),  The  corresponding  values  (a,R)  are  then 
determined  by  means  of  the  transformations  (6.8).  The  points  (a,R) 
for  all  possible  values  of  c give  the  curve  of  neutral  stability 


in  the  (e»R)  plane. 

lividently  z(c)  and  u(c)  could  be  obtained  graphically  as  just 
described  by  drawing  graphs  of  l(L,(z,c)  end  (z,c)  against  z for 
each  value  of  c.  However,  for  the  purpose  of  the  numerical  calcu- 
lations it  is  more  convenient  and  more  accurate  to  solve  the  two 
real  relations  obtained  from  (6.1)  for  z and  u by  an  iteration 
procedure  with  the  aid  of  the  tabulated  values  of  5„(z)  and 
(see  Table  1),  as  described  in  [2].  The  maximum  value  of  c and  the 
corresponding  vrlues  of  z and  u can  also  be  found  accurately  by  an 
iter? tl on  method  (see  Appendix), 

Present  calculationr  are  carried  cut  for  the  boundary  layer  on 
a flat  plat®  at  a Mach  number  of  M-  = 1.6  and  ratio  of  wall  bo  free- 
stream  temperature  of  = I.O73,  This  temperature  ratio  is  the 

critical  ratio  below  which  the  boundary  layer  is  completely  stable 
with  respect  to  two-dimensional  disturbances  (see  Section  10).  The 
Prandtl  number  is  given  the  value  0.7^,  and  for  simplicity  of  cal- 
culation the  viscosity  coefficient  is  teken  to  be  directly  propor- 
tional to  the  absolute  temperature,  with  the  constant  of  propor- 
tionality chosen  so  as  to  ^Ive  the  correct  result  at  the  wall. 

This  procedure  for  the  vlsooaity  coefficient  is  known  to  be  a good 
approximation  in  the  case  of  the  boundary  layer  itself  when  the 
Mach  number  is  not  too  large.  Although  it  probably  is  not  as  accurate 
in  the  stability  calculations,  estimates  shows  that  it  does  not 
algnlfiaantly  alter  the  general  trend  of  the  results <,  A detailed  dis- 
cussion of  the  calculations  is  given  in  the  Appendix. 


In  Tajle  3 a,a^,aQ,  R,R^,Ro  ape  ta'ouletcd  ts  functions  of  c 

A ^ V X w ^ 

for  several  values  of  In  th&  ranp,e  0 < < 1,6,  vjhsre  a_ 

and  R are  wave-nuriber  and  Reynolds  number  based  on  a reference 
length  Xj^{the  dlstrnoe  from  the  leadlnj-j  edge  of  the  plate)  end 
ttg  and  Rq  are  those  based  on  the  momentum  thickness  as  reference 
length. 

In  Figure  3 neutral  curves  in  the  (a^.R^)  plane  are  shown  for 
the  three  cases  = 1,4>  1.0, 0.6  corresponding  to  direction  angles 
of  29.0®,  51.3^,  68. 0'^  respectlvc;ly . 

The  approximation  is  made  that  the  minimsam  critical 
Reynolds  number  is  very  close  to  the  value  of  R corresponding  to 
the  maxlinum  value  of  c.  This  uos  found  to  be  true  in  previous  In- 
vestigatiens  for  an  incompressible  fluid  [10 ] and  for  e com- 
pressible fluid  [?j.  The  neutral  curves  1n  Figure  3 show  that  in 
the  present  case  it  is  also  an  excellent  approximation,  with  only 
a very  slight  error.  The  minimum  critical  Reynolds  number  approxi- 
mated in  this  way  is  plotted  against  the  direction  angle  of  the 
oblique  waves  in  Figure  4 and  tabulated  in  Table  4*  might  have 
been  expected,  becomes  infinite  at  direction  angles  of  0°  and 

90  , i’he  minimum  value  of  Rx_  occurs  at  M-i  = 1.0  30  (corresponding 

C •*'  { '\  ')  A 

wC  s clijTOC tioii  oi*  is  i^iVeri  lsv  Hxq  ' ^ “ 1 • 6i4-  ^ —0  * 

This  value  Rx  = i-,64  x 10^,  therefore,  is  the  minimuin  critical 

c 

Reynolds  niamber  for  all  three-dimensions  1 disturbances  in  this 
particular  case. 

9.  The  dependence  of  the  minimum  criticel  Reynolds  number  on  the 
direction  of  propagation  of  the  waves 
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For  those  cases  in  which  the  Eppreximations 
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(8.16b)  are  sufficiently  accurate,  it  is  p-Cnsible  to  derive  fairly 
simple  analytical  results  for  the  directional  dependence  of  the 
minimum  critical  Reynolds  number  (that  is,  dependence  on  the  para- 
meter ) , By  solving  (8.16a)  for  a we  have 
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Then  the  minimuji  vxlue  of  R on  the  neutral  curve  in  the  (a,R) 
plcne  is  ^l^/en^  with  the  aid  of  (6.6)  snd  (8,2a),  by 


(9.3)  Rt 


I t Ji  H,  (Xf, 
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whore  c la  approximated  by  the  naxlinu;n  permissible  value  of  c 
described  previously.  Cncathe  particular  boundary-layer  fbw  Is 
specified  (i.e.  the  functions  w(y)  and  T(y)).thls  value  of  c 
and  therefore  z(c),u(c)  and  all  the  other  functions  of  c in  (9.3) 
are  det^L,^ralned.  For  small  X,  previous  discussions  show  that  c la 
-iven  approximately  by  v^( c )r:o . 5‘Ou * end  also  that 
z(c)  ~ u(c)  2:  I.I4.8  (see  Figure  1). 

U‘e  now  have 


where  K depends  on  the  given  boundary-layer  flow  but  is  independent 
of  In  order  to  detarralne  the  dependence  on  of  the  other  quanti- 
ties entering  Into  R^,  we  note  that 


(95)  H, 
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la  Independent  of  Mt,  and  that 
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Thon,  with  the  al<]  of  the  definition  for  k,  (9.^4-)  can  be  written  as 


(9.7)  tv  - 


A tBM, 
n.tji-FiK'-c)^ 


where 


(9.8a)  A ^ K fO-<^>^j  + r~J 

(•  .J  (i-crj 

(9.0b)  B - K jO  ’^)  ~ i 

The  quantities  A,B,o  are  independent  of  and  are  determined  from 
the  distributions  w(y)  and  T(y)  of  velocity  and  temperature^,  respective" 
ly,  in  the  boundary  layer. 

From  ( 9 ..  7 ) wo  find 

(9  9,  d Rc  . 
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la  the  condition  that  ell  throe-dlmenaional  dlsturbsncoa  have  a 
hlKhsy  mlnirauia  critical  Reynolds  number  than  tvio-diiuens j.cnel  ouQS. 


Thus,  if  it  can  bo  shown  that  the  parameter  C 


- Ml  /2( 


i-c)  = 


...  .. 


for 


a epecified  boundary -layer  flow  is  less  then  unity,  then  three- 
dimcnslonsl  dlstui-bances  sre  tnure  stable  and  are  not  very  Iriipori'.ant 


in  the  stability  problem. 

In  Table  5 Figure  5 given  the  results  of  calculstit.'ns 
of  C S3  a function  of  for  an  Insulated  i'i.at  pl^te  from  = 0 
to  Mt  = 2.4,  and  of  C as  a function  o;(  Tj^/T^  at  o Mach  nuraber  of 
Ml  = 0.7  from  ^ 0.?0  to  T„/T2  1.23'.  For  the  tiisrulatod 

pl.ato,  it  is  cleai>  that  three-dimensional  diatui'baness  ar©  more 
stable  than  two  up  to  a Mach  number  of  about  i«8  they  ba- 

cc2i«  aiightly  more  unatoble.  The  f ict  that  fl  decreases  and 
again  becomes  loss  ths.n  unity  %.'noa  It  increases  beyond  lo8 
may  not  be  of  muen  signi..i  loance  oecause  of  the  large  errors  at 
such  Mach  numbers  due  to  the  approxiIjnatiouO  used  In  the  deductions. 

10.  The  Influence  of  surface  cooling  on  boundary-layer  stability 

The  problem  of  boundary-layer  stability  at  supersonic  speeds  is 
difforent  from  that  at  subsonic  speeeja  in  one  important  regnect-there 
is  always  a class  of  oblique  waves  {•lliose  with  1 < ^ M-,  ) for  which 

the  values  of  c are  bounded  away  from  zero  (i.e.  c > 1 - 1/Mq^  >0). 
This  leads  to  interesting  consequences  in  the  case  of  surface  cooling. 


Let  U3  first  consider  strictly  tv;o-dimensional  disturbances 
(for  which  since  the  situation  is  then  somewhat  simpler. 

'ihe  irilriimum  value  of  c on  the  neutral  curve  is  1-1/M^,  and  the 
intxiraura  value  is  ^;iven  by  v^(c)=  t.  (c),  that  is,  it  is  the  value 

O Ijfl 

of  c corresponding;  to  the  intersection  point  of  the  (;reph  of  v^(c) 

with  that  of  Ipi  (c)  (see  i’if;ure  2).  For  sufficiently  small  T /'f.,  . 

m " 

v^(c)  > 0 for  c = 1-l/M]^  and  Vq(c)  la  monotonically  increasini;  for 
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the  Vjj  curve  higher  In  the  (CsV^)  .plane.  This  behavicur  of  the  v^^ 
curve  is  typical  for  Mach  numbers  up  to  For  Mech  num- 

bers greater  than  = 7-46  no  amount  of  cooling  is  sufficient  to 
raise  the  curve  sufficiently  to  make  v^(o)  positive  at  c = 1-l/M^, 
We  see,  therefore,  that  in  this  range  of  iiach  numbers  for 
T /T,  a. nail  enough. the  range  oi  values  oT  c on  the  neutral  curve 
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an  be  made  as  c'loae  to  as  desired 


by  maKlng  smaller.  Thus,  there  is  a critical  toiTipereturs 

ratio  (T^T^)^  given  by  v^(l~l/Mj^)  = ) for  which  the 

whole  neutral  curve  disappears.  This  can  be  seen  by  solving  (8.16a) 
for  a and  noting  that  a-*>  0 os  c-f  1 . Therefore,  when 

(TyTi)c,  al  1 v^eluea  of  c on  the  neutral  curve  approach  l-l/H^j 

/V 

so  that  all  values  of  a~>  0 and  all  values  of  R->°^. 

It  is  to  be  noticv^d  that  in  this  limiting  situation  relations 
(6.16a)  and  (8.16b)  becomes  exact  since  the  errors. are  0(a*^),  so 
that  the  relttlon 
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the  value  of  the  Prandtl  nuxnber  end  the  particular  viscosity  lew  used. 
The  value  M-  = 7*^!-6  corresponds  to  a I’randtT  number  of  b.7‘^  and  a 
linear  viscosity-temperature  relationship. 


is  oxeet. 


for  t.ie  critical  tci’ipera  ture  retio  ( T /T- ) ^ 

Vheu  evidently  v^(c)  > c > 1-  1/M^ 

so  that  no  solution  is  possible.  Therefoi-e  for  ^ V 

the  boundary  layer  is  completely  stabilized  with  respect  to  tvio- 
dimenslonsl  disturbances. 

Calculrtlons  (doscribod  in  detail  in  t ie  Apijendix)  have  been 
carried  out  to  determine  the  variftion  of  (Ty/Tj^)^  with  for  the 
boundary  layer  on  a flat  plate  ■with.  Prendtl  number  ecjuel  to  t'  = 75 
end  viscosity  coefficient  directly  proportional  to  absolute  tempera- 

^ . - — - ^ > |l  t*.  ^ ^ ^ t ^ ^ ^ ^ ^ 4 «.»  W ^ ^ ^ ^ ^ M « «ta»  Vr  r.\  ^ ^ ^ ^ 
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cal  te.riperature  ratio  against  Mach  nu;.'iber  for  two-dimensional  dis- 
turbances is  shown  in  figure  6.  Since  the  critical  temperatu-’e  ratio 
approaches  zero  at  Mi=  7.46,  no  amount  of  cooling  la  sufficient  to 
completely  stabilize  the  boundiyy  layer  with  respect  to  two-dimen- 
sional disturbances  when  M-  > 7.I4.C.  The  present  mumerical  results 
do  not  differ  eopreciably  from  those  of  Vtn  Driest  [11],  even  thoup;h 
corrected  viscous  solutions  (which  enter  into  the  detev*mlnetlon  of 
X(c)  and  thus  of  fijj,(c|)are  used.  The  explanation  is  prooebly  the  fact 
that  v^lc)  is  so  extresasly  sensitive  to  variations  in  Tw/'^X  ssakLl 

errors  in  in  (10, i)  do  not  affect  the  solution  very  rnuch.^’ 

The  conclusions  depend,  of  course,  very  sti*ongiy  on  the  value 


oi*  niint'SX^  s.nci 
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Irtlonship  used  (see  [llj).  However-  for  the  same  re;  son  as  given 
above,  it  is  not  likely  that  the  use  of  the  present  viscous  solutions 
would  significantly  alter  the  results  for  these  other  ersos  obteinod 
in  [111. 


*Van  Driest* 3 computations, in  addition  to  the  use  of  the  Inaccurate 
X(c)  of  reference  [Ij,  also  inv'-olve  what  is  essentially  an  approxima- 
tion for  fim(c)  in  terms  of  X(c).  It  is  interesting  to  observe  that 
if  the  exact  procedure  of  calcu  Xs.  b i.  CT3  is  used,  in  combination  with 
the  inaccurate  value  of  X(c)  ’.iven  by  reference  [l], several  critical 
cooling  retivS  will  ue  found  for  each  value  of  c=l-  This  has 

been  pointed  out  by  Professor  Martin  Bloom.  However,  wluh  the  pre- 
sent function  X(c)  wnose  Influence  is  small,  such  difficulties  do 
not  arise. 


I-iow  consider  oblique  w£vea  with  direction  angle  arc  cos 

where  1<M^<M^=  Just  as  in  the  case  of  strictly  two-dimens ional 
disturbances,  there  is  again  the  possibility  of  completely  stabi- 
lizing the  boundary  layer  with  respect  to  oblique  waves  proi.-a  gating 
In  this  direction  for  small  enough  T^/T^.  The  arguments  carry 
through  very  much  as  before  except  for  one  co.nplicatlon  due  to 
the  fact  that  v;hsn  the  Jisch  number-  Is  Isr-ge  the  curve  Vq(c) 
against  c may  have  an  inflection  point  and  thus  a maximum  and  mini- 
mum In  the  inucrval  0 --  u ^ 1 (see  Figure  2(b)).  Thus  even  when 
T^T-,  ia»fflBll  enough  for  the  relation 


(10.2) 


v.d 


to  hold,  the  corresponding  v (c)  may  not  be  monotonlcally  increas- 
ing  for  c > 1 - 1/Hjf  when  M,  la  small  enough'' . 'n'‘oro  could  still 
be  values  of  c > 1 - 1/M^^  for  which  v (c)  < (c)  (see  Figure 

2(b)),  so  that  a neutral  curve  might  still  exist.  The  critical 
temperatui-'e  ratio  below  which  the  boundary  layer  is  completely 
stable  wich  respect  to  oblique  waves  with  directlori  angle  arc  cos 
M-j^/Ml  is  therefore  the  largest  T^/Tj  for  which 
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Tnis  value  of  T,^/T^  Is  evidently  the  one  "or  which  tijs  cerreapondin.g 
v^(c)  curve  Is  just  tangent  near  its  mlnlmu  i point  to  the  i}lijjj(c) 
curve  (see  Figure  2(b)). 

Obviously  if  W,  < 1 the  minlmuiri  value  of  c on  :he  neutral  curve 
is  just  c=  0.  Therefore,  it  is  not  possible  by  any  amount  of  cod- 
ing, to  completely  stabilize  the  boundary  layer  with  respect  to 
oblique  waves  with  direction  angles  greater  than  arc  cos  l/.M^.  The 
limiting  directions  are  thus  the  normals  to  the  Mach  linos  In  a 
plane  parallel  to  the  solid  surface. 

For  the  boundary  layer  on  e flat  plete  at  s Mach  nuirioer 
= 4 with  Prandtl  number  equal  to  0,75  end  a linear  viscosity 
t6'-iper£ ture  relation,  the  critical  temperature  ratio  T^/T-j^  has 
been  calculated  as  a function  of  direction  angle  arc  cos 
and  plotted  as  a full-line  curve  in  Fl^'ure  7.  The  cal  r.m  nti  ons 


Tnis  clfficuloy  does  not  occur  when  Mt  Is  as  large  as  Mi, at  ler.st 
for  the  boundary  layer  on  a flat  plece."^ 


>» 
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ar«  carried  out  by  first  firidi.ig  T^/T^  ss  a function  of  from 
(10./3).  Tile  resulting;  values  of  T /T,  when  plotted  sfrainst  arc  cos 

s curve  with  a maximum  and  e.  minimum  for  in  the  ranj^e 
0 < < M^.  The  proper  curve  for  the  critical  temperature  ratio 

is  the  obtained  by  drawing  a horizontal  tangent  to  this  curve  at 
its  mlnlmutn  point  (thus  cutting  off  the  dashed  part  around  the 
maximum  point  In  Figure  /)  . The  calculation  procedure  is  de- 
scribed more  fully  in  the  Appendix,  and  the  numerical  results 
are  listed  in  Table  7» 

Figure  7 shows  that  for  a temperature  ration  T^/T^  = I.70O 
all  two-dimensional  disturbances  ere  completely  stabilized,  and 
for  a sOTiewhat  smaller  ratio  = 1-474  the  three-dimensional 

disturbances  corresponding  to  obliaue  waves  with  direction  angles 
from  0^  to  74-3°  sre  completely  stabilized.  The  critical  tempera- 
ture ratio  decreases  to  zero  very  rapidly,  however,  when  the 

direction  angle  increases  from  74-3°  to  arc  cos  1/M^  = 75- 5° • 

Vvhen  the  tempei*ature  ratio  T^/T-j^  is  as  low  as  1.474»  there- 
fore, all  oblique  waves  are  stabilized  exceot  those  with  dlrectioii 
angles  greater  than  74 -3^-  iistlmates  of  the  minimum  critical  Reynolds 
number  of  the  letter  (all  of  which  have  very  sinall  values  of  th« 
wave-spssd  c)  show  that  10°  or  Rxq'^  10^  * Since  such  large 
values  of  the  Reynolds  number  are  far  beyond  the  range  of  practi- 
cal importance,  it  follows  that  if  the  wall  temperature  ratio  can 
be  reduced  to  as  low  a value  as  1.4?4  the  boundary  layer  will  be 
very  highly  stabilized  from  the  practical  point  of  view.  Such 
cooling  is  fairly  extrema,  however.  The  situation  for  larger 
values  of  T^T^^  is  not  clear,  except  that  it  seems  likely  the 
boundary  layer  is  quite  sensitive  to 


This  is  equivalent  to  the  use  of  the  condition  mentioned  in 
connection  with  relation  (10.3)- 


changes  In  a small  increase  above  = 1.4?4  niight 

make  the  boundary  layer  quite  unstable. 

Kor  each  value  of  up  to  = 7*48  i-3  probe blo  that  tlie 

naturo  of  the  oblique  waves  Is  about  the  same  as  at  = I4.:  that 
la,  for  a sufficiently  small  temperature  ratio  it  la  possi- 

ble to  completely  stabilize  oblique  waves  with  direction  angles 
less  then  arc  cos  1/M-,  but  not  those  with  lar.?jer  direction  angles. 

X 

It  Is  clear  from  Figure  6,  however,  that  such  temperature  ratios 
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At  the  present  time,  very  iittie  more  can  be  said  about  the 
general  problem  of  boundary-layer  stability  at  such  large  Mach 
numbers.  The  complicated  shape  of  the  v^\c)  curve  in  many  cases 
of  Interest  makes  the  problem  difficult,  as  we  have  seen,  and  also 
the  present  methods  of  numerical  calculation  do  not  appear  to  be 
adequate  except  for  the  unimportant  class  of  problems  involving 
extreme  cooling  (when  a Is  very  small). 

For  lower  Mach  numbers  (up  to  about  2),  the  picture  is 

exuaror*  iuti  uicax  rt^^iuxts  oooi^xuu  c?i'o  oncaa 

factory  approximations,  provide  considerable  useful  informetion. 
The  numerical  calculations  indicate  that  the  parameter  G increases 
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plate  at  a given  Mach  number  is  evidently  to  make  three-dimensional 
disturbances  relatively  more  unstable,  so  that  with  sufficient 
cooling  some  of  them  can  actually  become  more  unstable  than  two- 
dimensional  disturbances.  Since  C for  an  Insulated  plate  increases 
from  C=0atM2^=0toa  value  close  to  unity  when  ~ 2,  It 
aouears  that  the  amount  of  cooling  required  to  make  any  of  the  three* 


dlmonslonal  dlsturbeniies  more  unstabls  than 


)“dimonsl 


extreme  at  low  subsonic  Mach  numbers,  but  gradually  becomes  less  end 
less  extreme  as  M^  increases  until  at  a Mach  number  of  the  order  of 
two  very  iittie  is  required.  In  particular,  at  M-j^  = 0.7,  the  wall 
temperature  ratio  must  be  reduced  below  = 0«58 

(see  Figure  p). 

It  is  known  that  at  subsonic  Mach  numbers  cooling  of  tiie 
solid  surface  is  very  effective  in  stabilizing  the  bomdary  layer, 
if  only  two-diraenaional  disturbances  are  considered  (sec  i2j).  The 
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valid  even  if  three-dimensional  disturbances  are  considered,  since 
ths  iattar  usually  have  a higher  minimum  orltiual  R_’'ynolda  number 


than  two-di.menalonel  dlaturbancea  (except  undex^  conditions  oi 
extreme  cooliru;). 

Even  at  the  moderate  supersonic  Mach  numbers  when  three-di- 
mensional disturbances  ber^.ln  to  play  the?  leading  role,  cooling 
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stabilization  at  supersonic  Mach  numbers  as  predicsd  in  [2}  is 
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Show  that  for  an  insulated  flat  plf  te  at  = 1.6,  Hx^  is  of  the 
order  of  10^,  whereas  the  celcuistcd  results  show  that  for  a 
cooled  plate  at  the  sarae  Mi£  ch  number  = 1-6^  with  a tempera- 
ture ratio  ~ i»073  which  is  Just  small  enough  for  two- 


dimensional  disturbances  to  be  completely  stabili 
the  order  of  10^. 
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il.  Conclusions 

(11*1)  An  order  of  m'xgnitucs  analysis  of  the  coaiplote  llue5;c5.zod 
equations  for  a thi’ce -dimensional  disturbance  superimposed  on  a 
two-dimensional  boundary  layer  shows  that  they  can  be  reduced 
to  much  simpler  equations  in  the  first  approximation.  These 
simijlified  «qr«t5_ona  ere  still  valid  at  high  f“se-stream  Mach 
numbers  If  correct  velocity  and  temperature  profiles  for  the  boundary 
layer  in  such  cases  are  used.  At  very  high  Mach  numbers,  however, 
the  accuracy  of  a first  approximation  may  not  be  very  great,  so  that 
higher  order  ajjproximatioria  and  therefore  more  complicated  equations 
with  mors  terms  retained  may  have  to  be  considered, 

(11,2)  The  present  equations  for  a three-dimensional  disturbance 
-transform  to  equations  which  have  the  same  form  mathematically  as 


those  for  a two-dimonslonal  disturbance,  but  which  do  not,  how- 
ever, corres'^ond  to  a roal  two  dlmeuaionel  disuuraance.  Althongh 
no  direct  conclusion  regarding  the  importance  of  three-dimen- 
sional disturbances  is  now  possible,  the  use  of  these  transformed 
two-diiaenslontil  equations  lasds  to  a considerable  raathemetical 
simplification  of  the  theory. 

(ii.3)  Corrscted  viscous  solutlcns  have  been  obtained  which  lead 
to  a condition  for  the  charactei'latic  values  of  very  much  the 
same  form  as  that  in  ths  previous  theory  of  Lees  and  Lin.  Detailed 
mathematical  discussions  will  be  given  in  another  report.  The 
present  formulation  of  the  theory  porm'ts  more  accurate  numerical 
calculations  at  supersonic  Mach  numbers,  when  the  wave-speed 
c is  often  large.  It  reduces  to  the  Lees-Lln  formulation  for  small 
values  of  c,  so  that  the  validity  of  their  results  in  such  cases 
is  verified. 

(n  Ji)  Consideration  of  the  present  analytical  and  numerical  re- 
sults sug-^ests  that  probably  three-dimensional  disturbances  ere  of 
little  importance  at  subsonic  free-stream  Mach  numbers,  since  they 
usually  have  a higher  minimum  critical  Reynolds  number  than  two- 
dimensional  ones,  except  possibly  uxider  conditions  of  extreme 
surface  cooilngo  However,  further  numerical  calculations  for 
oblique  waves  in  this  range  of  Mach  numbers  would  be  useful. 

(11.5)  Mach  number  increases,  three-dimensional  disturbance 
enter  the  picture  under  oonditious  which  become  less  and  less 
extreme,  until  finally  for  Mx  between  one  end  two  they  begin  to 
play  the  leading  rols  in  many  problems  of  practical  interest. 

(11.6)  At  supersonic  free-stream  Mach  numbers  the  boundary  layer 
can  never  be  completely  stabilized  with  respect  to  all  three-di- 
mensional disturbances.  Although  oblique  waves  whose  directions 
of  propagation  are  at  an  angle  leas  than  arc  cos  l/M^  with  the 
free-stream  direction  can  be  completely  stabilized  by  sufficient 
cooling  when  is  not  too  large,  those  whose  directions  of 
propagation  are  at  a greater  angle  cannot. 
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(11.7)  The  graph  of  critical  temperature  ratio  (belou  which 
strictly  two-dimensional  disturbances  are  completely  stabilized) 
versus  Mach  number  (Figure  6)  is  not  significantly  different 
from  that  found  by  Van  Drleet  [11 3,  In  spi.to  of  the  corrections 
to  the  viscous  solutions.  Although  these  corrections  ere  very 
important  for  the  accurate  determination  of  the  curve  of  neu- 
tral stsbliity  at  supe'-'scnic  Mach  lamabex’Sj  they  are  unimportant 
for  the  determination  of  the  critical  temperature  ratios. 

(Ilb8)  For  i'lach  r.uribers  up  to  abotit  t-.T,-.;..  cooling  of  the  solid 


surface  is  very  effective  in  stabilizing  the  boundary  layer,,  At 
higher  Mach  nuiubaro  it  appears  that  the  arjiount  of  cooling 
necessary  to  ralsci  tho  minimum  critical  Keynoldr,  number  signifi- 
cantly becomes  pr-acticslly  prohibitive. 

(11.9)  Very  little  more  can  be  saia  about  the  problem  of  boundary- 
layer  stability  at  Mach  numbe?s  greater  then  two,  since  the  present 

of  numerical  caiculaticn  for  tho  curve  of  neutral  stability 
(and  thus  tho  minimun  criticfBl  Reynolds  number)  do  not  appear  to 
be  adequate  in  s,:cb.  cases.  The  main  difficulty  is  in  evaluating  the 
inviscid  solitbions  of  the  disturbance  equations  when  e and  a are  both 
large.  This  pi'oblern  is  being  studied  further. 

(11.10)  A more  serious  difficulty  then  the  preceding  Is  tliat  ac  very 
high  Mach  numbers  the  present  asymptotic  theory,,  which  is  essentially 
a first-order  theory,  may  In  itself  by  Inaccurate  (as  pointed  out  in 
(ll.l))o  More  inf oriiiation  about  this  matter  will  be  provldod  by 
numerical  resul bs (including  estimates  of  certain  parameters ) for  high 
Mach  numbers,  vaaich  will  be  obtained  as  noon  as  the  ccruputatlonsl 

1 1-  ^ r>  e yr»or>  -f*  -iio  (11  d \ <>.  /^1  {-9 
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(11,11)  It  is  believed  that  tho  approiiimets  method  of  calculation  of 
the  neutral  curve  described  In  Section  8 Is  adequate  up  to  e Mach  number 
of  the  order  of  two*  eapeclelly  In  cases  Involving  cooling,  when  r and  a 
are  both  small.  Addltioncl  nautr^^  1 curvos  for  the  oblique  waves  In  this 
range  of  Mach  numbers  or  perhaps  further  c«l0ul*tiecui<Ctiak  typ«  dMKnrib^di 
In  Sections  9 and  10,  woulo  be  of  considerable  value.  It  might  also 
be  of  Interest  to  determine  In  a few  typical  cases  the  effect  of  using 
a more  accurate  approximation  than  f8.l6  a),  by  keeping  terms  of  higher 
order  In  a,  as  Lees  has  d >ns  1 2], 


r or  prise  L- 1 c ti  j,  pur[>otst!i:> , the  transit  lOn  IvSj^nolds  nuiiiber  Is  muen 
more  Important  than  the  minimum  critical  Reynolds  number.  Ihus, 
Information  concerning  the  rate  of  amplification  of  the  boundary-layer 
disturbances  after  they  enter  the  unstable  region  would  be  extremely 
useful.  Further  investigations  are  being  directed  tovjerds  developing 
the  amplification  theory  for  three-dimensional  diaturbances , vihich  mpy 


be  different  from  that  for  strictly  two-dimensional 

disturbances. 


Appendix*  Procedure  for  the  numerical  cslculatlone 

(a)  The  traneforraatlon  to  w a£  Independent  variable 

Moat  cf  tho  stability  calculations  can  be  carried  out 
moat  conveniently  with  w as  Independent  variable.  This  formu- 
lation has  also  been  found  to  be  very  useful  for  the  boundary- 
layer  calculations  (see  [12j  and  I13])* 

First.  Introduce  new  non-dimensional  variables  as  fellows: 


it\ 


t 


r> 

Y 


»IV 


/t 


Where  'f  « nw>  and  w.T.p.  n are  as  given  in  the  List  of  Symbols. 

f / 

Then  relations  (6.2a)  and  (8.2b)  can  be  transformed  to 


(2a) 


- - U) 


<X  M = 


= 


^ / bl.  r * "T  \ 
- IT) 


r(t 


j 


-L  r/c)  - 1 
2 ' 


Where 


(3) 


1(c)  = ^ 


Also,  relation  (8.1b b)  for  v^(c)  becanes 

o 


*The  present  symbol  6 la  rot  to  be  confused  with  the  same  symbol 
used  slsewhara  in  this  r^nort 


-1^0- 


where  ©,1J>  are  considered  aa  functions  of  w and  ,s' 

ore  their  derivatives  with  respect  to  w. 


The  non-dimensional  vlscosltv  coefficient  V la  related  to 

I 

6 most  accurately  by  the  Sutherland  viaccslty  law  (see  I13],  p«7) 

G’S'OS' 

7;;. 

O'S-QT 

The  non-dimensional  tempersture  6 Is  given  as  follows 
(see  [13]*  p»6)  when  the  specific  heats  C_  and  C_  are  constants: 


(5) 


^ = e 


% 


/ + 


e + 


(6)  $ - \ - (i-  = %(^) 

/<*>  * TV 

where  Y * ®p/^v  “ functions 

tabulated  on  piZy  of  [13J*  CrOcsO  [12]  hsa  shown  that  these  func- 
tions are  practically  Independent  of  the  particular  vlscoslcy- 
tssiperature  relationship  used. 

For  the  case  of  an  Insulated  wall,  the  ratio  of  wall  to 
free-atream  temperature  is  given  by 


(7)  T„  = J+  l/F  ^ 

/ \ 

\ ; 

p it  follows  from  (1)  that  ^ = 9o 


When 


P/i 


« constant 


rA 


The  function  s(w)  = gfl.('w  )/g^.^(0 ) la  detemalned  fi-om  th© 
function  g|j(w)  «hlch  satisfies  the  following  differential  equation 
(see  1133#  p.6); 

- o 

with  boundscy  conditions  g|.{0)  = O5,  g^(l)=0<,  If  p ^ = ^w/S# 

• then 


^*9'  ... 

( 9 } “ • T ^ ~ u 

K K 

where  - p„  Therefore,  g (w)/K  = g„(w),  the  Blasiua  shear- 
stress  function  (1.. o.  + 2w  ~ 0 v^ith  g^(C)  = Op  g^^(i)  = 0)c 

This  function  is  tabulated  in  [133,  p»25o 

For  p yU  = Py  therefore,,  s sntiafie.s  the  following  equnbion- 


(la) 


S s"  i-  -4r~  ^ -O 


with  boundary  conditions  s»  (0)  = 0,  s(l)  = 0,  since  s g4,(w)/g^(0) 
s=  g^(wj^g^(0)«  A coTiparison  of  this  Blasius  distribution  for  a with 
the  distributions  obtained  in  [I3]  when  the  accurate  Sutherland 
viscosity  law  is  usod  shows  that  in  many  coses  it  is  actually  quite 
a close  approximation  to  the  correct  one. 

The  function  6 satisfies  the  following  differential  equation 
(seo  1133,  p.6)? 

/a  t . \ 


(11)  I a -t  tr(  j jS  0-(t)  e >S’' - O 


with  the  boundary  conditions  9(0)  ~ 1„  6{1)  - 1/T 


I 
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The  Amotions  s(w),  «»(w),  GjCw),  ftj{v),  6jj(w), 
are  tabulsted  In  Table  8.  The  values  for  and  6^^  have  been  re- 
calculated by  the  method  described  in  [13]»  and  are  some^at 
more  accurate  than  the  values  tabulated  in  [13 ]»  p.27«  The  Prandtl 
riuniber  ^ Is  taken  to  be  Of/5* 

Per  small  values  of  c«  1(c)  and  X(c)  may  be  obtained  most 
conveniently  from  power  series  in  c.  With  the  aid  of  equations 
(2b),  (3)  and  (10)  and  the  fact  that  g*(0)  » 0.66411,  these 
series  are  detemlncd  as  fellows  (for  p yu  = P,,yUy)t 


(12)  Ifc) 


£ 

3 


i7T67S\jiO'U^li)*-\ 


(13) 


Me)  - 


3/r  j 


3 I 
C H- 


iii-r_L_Tc 


=:  0^n5'l6j  C ^ O-O^'fOOS  C + 0(c  ) 


The  functions  I(o)  and  X(c)  are  tabulated  in  Table  8„  the 
values  for  small  c being  obtained  by  use  of  (12)  and  (13)  and  those 
for  large  o by  numerical  integration  of  (3)  and  use  of  (2b). 

In  the  calculations  the  quantities  (S/x-^y/^  •‘Od  T/ 
are  reqiilred.  In  the  Croeco  method  the  velocity  distribution  is 
given  by 


z C 


i 

j 

O 


rx 

• JT 


(see  [13j,  p»20)  sc  that 


(15) 


2 


M/r 


L \JJl  = 

«■!  J % 


Where  w-,  la  the  velue  choaea  to  define  the  edge  of  the  boundary 
layer..  In  the  preaent  ealoulatlons  the  value  of  is  taken  to 
bo  0*999«  Por  p yu  = equation  (15)  becomes 


<161  i ifT' 


- 9 


( ^ T,/j  r.j 


= 2 3-00  - 2 OT90(i-  j:  U o <.f<r»r>'-<)W,’‘;~  ^ 

? i aV'  ^ 

Pra.n  equation  (l^-O)  p.  17  of  (133  we  see  that 


/ ^ 


(17) 


1 r 

/V  — — — ( ^ 

/5-.'U,V 
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-L 

_ t 
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* i/^; 
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/X, 

I S I T"  \ 


Therefor# 


nflA 


Tv  . 


The  value  of  yu^  la  determined  from  the  Sutherland  viscosity  law 


‘■kh- 


( > 


(»••  [133  p.  7): 


(19) 


y*. 
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(b ) The  tgensf owne tlon  to  m Independent  variable  In  the  outer 
part  of  ths  bound^&sr?  layer. 

In  the  evalxiatlon  of  certain  Integrals  that  are  required  In 
the  calculations  It  is  more  convenient  to  use  some  Independent 

VCIX'XCLUXC  WuAoF  UXlOil  ■■  XIX  VXXO  VUUOX*  ^OX*U  WX  UX4.0  UWUlXtXOA' jr  Vpo  AX4 

the  present  case  when  p yu  = ^v/V*  norrnal»dlstancs  pararrieter 
^ of  the  standard  Blaslus  solution  for  the  non-dimensional  stream 
function  f(i^  ) Is  suitable,  where 


(20>  MT  ^ I 4*(^)  . 

2 If  ~ ' 


and  44"=0 

p St 


wi  th 


boundary  cendltlcais  f(0)  « f*(0)  « 0,  f*{<»)»2  (see  [ilfj,  p->  135)' 

— « 


For  tho  BJLuoxUo  case,  yu 
Pl»  i ■ ^U^/*P,*Xj^)^3C2  and 


“**  constant  * M*p  P*  = constant 


(21)  ^ ± f(s)(  V,/v*M.f£  ^ if  i.  ^ /(t) 


But  from  (!T) 


(22) 


T= 


— * — O y _»  — s 

»^w  /^w ' ►'I  ri 


^VkAVSA 


U3)  h(^> 


/ A 


The  functions  0j(/  ).  ^(  j )s  «»(<“  )*  9Jj( 
values  of  may  be  obtained  with  the  aid  of  the  following  asymptotic 
formula  for  ^ f * ( ^ ) s 


(21^) 


i iw  ■ 


oo 


where  m » 1»73»  n = 0,231  (see  [153*  p®‘37)« 

7^(w)  «uu  Sjj(w)  are  dei'ineci  on  pp,  T-Q  of 


TV.a  -- 

— ebvs«av»  vxv/iAO 


[13]  as  follows: 

I(m^) 

(25)  0r  ('•*^)  = f77r 

^ V’/ 


o^ferM  JT')  -Jh")] 


where 


fry  / — 

IM  = \ £^:i'  <<«r 


(27) 


VUisn  the  transfon’iation  to  f as izidepexMleat  variable  is  made, 
the  following  asymptotic  forraulas  may  be  obtained  with  the  use  of 
(2L)I 

00 

.f  f 

(26) 


Pxi 


. 10)  l/F  J ^ 

Dir(f-a) 


(29) 


(30) 
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!~<T  ^ 


-o-f/-  f j- 


e'(()  - [/»]  c 

^ i LO) 


OO 
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^ i -t 


— •/  '* 


<Ir'(t-  f) 


CO 


H^/n' 


|/i(2-(r) 


|4a-(r)\  < 'f(j-<r/  f,*-ff  J ? | 


Ji  ({-«■' 


z / 


' J \ 

«TT*  / * ^ “•  • 

vzu-;^;  J 


Mher®  cr=  0o?5*  KD  = 1.1070,  f"(0)  = 1.32622,  J(l)  = 0„5768 
and  0.9252o  Then  ©jj(  j ) and  ©!tC^  ) from  (25)„ 

Prom  the  preceding  formulas  0^(5“),  ®j(  ®jt(^)»  ®ii^  S ^ 

now  bo  calculated  as  functions  of  ^ , with  the  eld  of  tables  of 
the  exponential  function,  error  function,  and  exponential  Integral. 
The  results  of  such  calculations  are  given  in  Table  8. 

The  valu^  of  the  Integral  1(c)  as  deflriKd  by  equation  (j) 
when  0*1  la  obtalneu  by  changing  to  ^ as  variable  of  Integration 
in  the  outer  pert  of  the  boundary  layer  as  fo3.1ows: 


(32)  !(•)= 
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r 0‘SS0  -\r0-0Z2-f0‘0&^  r 


Th«  fivalnniil-Tn  o?  intr^^^rsls  CCCV.2*irig 


/ T 4»  N 4 >*  vt  T 

y •dm ’*>»  f 


car'jpiod  out  most  tssily  by  ohenglng  ho  J.  neap  the  outer  edgft  ot 
the  boundary  layer: 

S(0'W) 


51AM. 
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ijjiMjr  ^ 


= \ 


(34) 
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\ ■ 0'^U-* 
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*• ' Otf'yS^ 
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ds  ^ Si^TIV  - 3^00 
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lliase  Integrals  are  evaluated  bv  numerical  integration  with  the 
use  of  the  functions  tabulated  in  Table  8« 

(c ) Calcula  tion  the  Integrals  H^(c ) and  Ki^(c ) 

The  integral  Hj^(c)  defined  in  (8o7a)  may  bo  expressed,  after 
several  transformations,  as  followst 

s f'ri  = ^ fr  ^(o)^-cfdt 

rr'  r.T.^j— 


\ t • M \ ^ / 


^{tL\  (Ln-^  U • io  330  -8-S^O  C H-i'OO 
\Ti^  / \ '^1  / \ / 


! 1 


o 


where  the  integrals  / w(  5 )d5  and  / w®(  ^ )d^  are  evaluated  by 


numerical  Integration^ 

Per  the  integral  k;,  defined  In  (8.7b)  by 

(37)  K., - n,  --  7k\_L_  ^ -Rr 

f J Cr^-c)»  ](^-c)^  ' 

o 

the  path  of  integration  goes  below  the  singular  point  y in  the 

c 

complex  y-plane.  Much  of  the  integration  can  be  carried  out  most 

conveniently  in  the  w-plane,  ao  that 

I ^ 

li/  ^A  ^ ^ ^ 


(38) 
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7W  J s(Mr-c)^ 
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M»-  i 


i j ,1 


Sfjur-c'i^ 
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where  0 < < w < w- ^ end  v-  = Oo999o 

GL  X X 

In  order  to  calculate  the  first  integral  expand  ©*/s  in  a 
power  series  about  w = c: 


(39)  = A 4 A + C t 0(MT^c)-t  ’ * 

where  A - B = (eVs)’,  c « |(eVs)Jp  T>  » i(e*/s)g,. 

Then 


(4o) 


^ -I  r_ 

\ ! - 1- 

\ LK-«)* 

M^cj  j ^ 
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Kio  i“ 

xne  value  of  is  evidently  given  by 

(41)  K,o  = -f?  -.  + B In  ^iZl 

c(Mi-c)  = 

The  other  Integral 

(42)  K„  - j [Ct  )+  • ' -]«U  = 

o 'b 

has  a regular  Integrand  for  ell  values  of  w.  The  upper  limit 
la  as  close  to  w^  as  the  integration  with  respect  to  w can  be  easily 
carried*  The  rapid  decreaBe  or  s towards  zero  for  large  w therefore 
limits  this  value*  In  the  present  calcu latioue  a value 
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w„  = 0.90  lf»  taken o 

Tha  coefficients  are  given  by 

p-(¥i 


U3)  i 


I C f i.%  iL^  - 2 f 'i'  + 


9 S 


■S  A 


-£l' 

2-  s J,. 


wl.'VwM  ciT  s «»ncA  6 at  w = o may  all  be 
obtained  in  terms  of  a,a‘,e,©»  at  w = c by  successive  differentiation 
of  the  differential  equations  (10)  and  (11),  as  follows; 


ujy 


\*t^/  ^ 


0"  r - (,.^)  e's 


- S S 2 


|.«a 


V &"  = --r '"«■■( (9 "+o-f »--()«, V )-h(i-trX9''s'i-  &' s") 


After  the  quantities  at  w ~ c are  determined,  the 

Integral  is  evaluated  by  numerical  integration..  The  value  of 
the  integrand  at  w = c is  given  by  G.  and  its  values  at  other 
points  are  obtained  by  direct  substitution  in 

A d 


iks:) 


5(U4t-c)^ 


In  general,,  it  is  found  to  be  possible  to  evaluate  the  integrand  In 
this  way  at  points  which  ere  close  enough  to  w = c to  permit  the 
integration  to  be  carried  out  with  sufficient  accuracy. 


The  final  integral 


(M>)  K,j,~  dur  = ^ (o)  — !21  „ olt  J/- 


\i^($yq 


»V  w/ 


X ^ \ 

5(-^; 


Is  easily  determined  by  numerical  lnte{0>ation 

The  functions  Kj^,(c),  ^igCc)  for  the  case  M,  = I06, 


» 1.073  •ro  listed  In  Table  8. 


(d)  The  determination  of  the  curre  of  neutral  stability 
The  relation  (8.1)  for  the  characteristic  values 

A.  V«  A ^ n-  i~i  #%  ^ A **l  A d.  X M 

uvf  WXAV  i#  isw  A'W^«awAwaA0 

<47*/ 

0 + 

(l  f 

fliTM 

% -r  • — 

= 

Vlth  V « 7q(s)  these  may  be  solved  for  s(c)  and  u(c)  by  a sudisessive 
approxiaations  procedure,  with  the  aid  of  the  values  of  ^^{z) 
and  Table  1,  as  fcllowas 


(l|Ba) 


<p.  - Cl±JLL 
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Slnce  X Is  usuelly  quite  small,  good  initial  approximations  are 
obtained  by  solving 

(49)  ^i(2-)  ~ yt/i(c)  , $^(z) 


From  the  discussion  of  Section  8,  It  follows  that  there  are 
usually  two  solutions  (s,u)  if  7^(0 } > 0 and  e is  not  too  large, 
and  that  there  is  a maximum  value  of  c for  which  tbese  two  solu- 
tions coincide  and  beyond  which  no  solutions  exist.  The  maximum 
value  of  e and  the  corresponding  value  of  z must  satisfy 


(50)  ^,(C)  - <p.(2.c)  ^ 


(i-i-  ;t)0! 


and 


(51)  « 
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d z 

^tieh  may  be  written  as 


(52)  3'  - 2.  A 

/ L'  • ''i'  - /«  y// 

The  aolutiou(z,o  ) of  (50)  and  (52)  may  be  obtained  by  the  success!^ 
approximations  scheme 
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with  tha  aid  of  the  functions  1^?^- r»5i  Table  1 and  the 
tabulated  values  of  Vq(c)  and  X(c)«  Initial  approximations  for 
c and  o are  given  by 


^ / 

(54)  “ 0 and  ^ r 


ie-  by  z=3»a2  and  M^(c)~0‘Seo^ 


With  z and  o determined,  u la  given  directly  by  (8«4a)« 

Whan  z and  u have  been  calculated  for  each  value  of  c,  a Is 
then  obtained  by  / solving  (8,l&a)»  The  quantity  u - L appearing 


4 A c«  as 


follows 


. wbjieh  follows  from  (£a)  and  (l8). 

1 

The  quantities  o and  R are  now  determined  by  means  of  the 

i 

transformations  (6,8),  and  tlienaj^,R^,aQ,Rg  are  calculated  by  means 
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For  the  case  = 1.6,  = 1.073»  t.he  parcaaetei*s^u^,  , 

V A-  •uiich  appear  In  the  calculations  are  given  by 

(59)  /iw-i-os<)s  , -6-3‘}6  ^ yk.  -i-isy 

r — 

(a)  Gclpiilation  of  the  pararostsr  C = Mj^y^vl-c )®-i- ^ 

The  values  of  the  parameter  C given  in  Table  5 have  bean  cal> 
culated  for  or  - 1 and  = 1,  for  the  sake  of  convenience,  alrico 
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then  many  of  tho  quantities  required  (sich  as  u-L  end  c)  can  be 
obtained  to  a sufficient  degree  of  approximation  from  the  data 
given  by  Lees  [2],  for  Mach  numbers  up  o = 1»3»  In  addition,  the 
boundary-layer  thickness  parameter  ~ I®  takaji  to  be  6.00  for 
tli«  plats  with  heat  transfer,  end  f'r  the  insulated  plate,  as  In 

[2j.  The  general  trend  of  the  results  would  be  expected  to  be  quite 
similar,  however,  for  values  of  c and  p 7loae  to  unity  and  for 
slightly  different  definitions  of  the  bovmt^ary-layer  thickness,  since 
the  various  calculated  quantities  are  not  v iry  sensitive  to  variations 
in  cr,  p yu  and  It  might  ba  mentionod  that  at  = 2.14.72  corre- 
sponding to  the  last  calculati^d  value  of  C 'or  the  insulated  plate, 
the  quantity  B/a  la  zsro  since  = (1-e'r  ir.  Ihia  case,  ao  that 
apparently  B/A,  which  is  positive  for  ibbbIIk' velues  of  ^ becomdB 
ns^atlve  for;*  larger  values. 

(f)  Determination  of  the  critical  iempcrsturs  x'Stlcs  f or  lii 3s t ion 

otythe  boundary  layer 

Prom  tho  considerations  of  Section  10  it  follows  that  the 
critical  temperature  ratio  in  any  particular  case  can  b e obtained 
by  solvizig 


(60)  yiC(o)  21  [j/-  ^ ^ z I- 

" Pf. 

for  the  ter.,  r.  ';..ture  ratio  T^.. 

I r^  \ 

For  given  c,  (53»)  with  0'“'  - ’3  solved  for  z and  then 

is  given  by  , v’aej  o is  defined  by  (8,4b). 

The  function  l|i  (c)  calculated  in  thiii  r.a.mer  la  tabulated  In  Table  6. 
* • 


Wh®n  p iU  = Pv/Hi  (60)  bsc 


ernes 


with  the  use  ol  reletlon  (ii.)  for  Now,  according  to  equation  (6), 
e Is  a linear  function  of  80  that  (61)  can  bo  reduced  to  a 

quadratic  equation  for  ao  follows: 


(62) 


T*“ 


Cii/’-  or^;^ 


= o 


(63)  ^ ’^(^"0^1  ^ 


The  quantities  6j(w),  6j(w),  6jj(w),  ©jj(w),  a(w),  s*(w)  are  all 

evslusted  at  w — c “ 1 ■=  1/mI.  The  coefficients  In  (6?)  nre  calou- 

„ % 

lated  with  the  asalatenoe  of  Table  8 and  the  equation  is  than 


easily  solved  for  T 


% O f ^ wsi««aa 


It  Is  to  be  noticed  that  the  v^alue  of  the  constant  C In  the 
viscosity— temperature  relation  y-^/  1 - C T / doss  not  enter 

into  the  formulas  (3)  end  (4.)  for  1(c)  and  Vq(c),  since  for  any  C, 
ip  = 6 for  this  linear  relation.  Therefore  ths  critical  tempera- 
ture ratios  do  not  dspend  on  C.  whirH  in  the  present  Investigation, 
it  may  be  recalled,  has  been  chosen  co  as  tc  give  the  correct  value 
for  the  viscosity  at  the  wall.  However,  in  the  calculations  for 
the  neutral  curve  the  value  of  C has  a small  influence. 
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It  should  also  be  mentioned  that  with  this  linear 
viscosity  relation,  the  functions  1(c),  X(c)  and  become 

universal  functions  of  c (for  a zero  pressure  gradient),  so  that 
the  values  of  these  functions  listed  in  Table  8 can  be  used  for 
all  values  of  ^ • 

For  more  complicated  situations  (Involving  the  accurate 
Sutherland  viscosity  laW,  pressure  gradient,  etc),  these  func- 
tions are  different  for  each  case.  However,  even  In  these  cases 
some  simplification  Is  possible.  In  particular,  from  equation 
(6«l4.b)  we  sea  that 

(64)  ^i(z,c)--  , '""f*  A-ACc)  j 

80  that 

^ 

where  is  the  maximum  of  Xj^{z»X)  for  a given  X.  Clearly 

^Ijq(X)  Is  a unlvei-sal  function  of  X,  which  rsay  be  calculated 

once  and  for  ail.  Thus  In  any  given  case  only  the  function  X{c) 

has  to  be  calculated,  and  the  value  of  1)fip,{c)  for  a particular 

c is  then  equal  to  value  of  X corresponding  to 

>y 

this  value  of  c.  A table  of  the  ftonctlcn  /vij^{Xj  coulcl  evidently 
be  obtained  by  just  re-tabulating  the  function  lTi^{c)  in  Table 
8 against  X, 


S~iX3  Xi 


}-g 

OL  i;y:riocl.s 


DijTienslonai 

quantities 

ilon-dim.Qnsionel  He  fere  nee 

quantities  quentlties 

Positional 

cw-tordinatea 

(1) 

X 

(2)  Xg 

y 

cJ 

(3)  X3 

z 

Time 

{hr) 

w 

1 

% 

Velocity 

coiTiponenta 

( ^ ) 11-  ~t  11  - 4.11  t 

' ^ ' ”1  ”1  ■ "J. 

- ( y ) ( j ) exp  f 1 ( ax+p  z-ac t ) 3 

(6)  Ug*  Ug+Ug 

V ( y ) +a4>  ( y ) exp  [ i ( ax+p z-ac  t ) 3 

(7)  u,  ^ uJ, 

h(y  )exp(  1(  ax+pa-ciot ) ] 

Density 

-Ox  ‘ft  —*  *•»*• 

v£)  p = p -J-  p 

P(y)"»-r(y )9xp[i(-7-x^^g^.^  act)l 

p. 

Pressure 

(9)  p'*^=  p*+  p** 

P ( y ) '1  ( y ) exrp  [ 1 ( ox+pz-ac  t ) j 

— ft 

p. 

Temperature 

( 10 ) T'- * T + T 

T(y)+i'*(y  )exp[i(ax+pz-act )] 

Viacoelty  0 oef f Ic  1 mta 

(11)  yu*.^*+yu’* 

‘ 'd’/i 

+ 6(y )exptl(ax-»-Bz-act )], 

f dT 

A 

C12)  X*-X*+>.*^ 

^(y)  ©(y)expIl(ax4pr-aot)) 

~e 

A 
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Dimensional 

quantities 

N o:a-  d ime  n a i 0 na  1 
quantities 

Reference 

Quentitles 

Thermal 

conductivity 

/i(y) 
o*(y  / 

(13)  k** 

^(  ^ ) 0(  y ) oxp  U ( Z' 

-act)] 

nfiVo—nUmbor  o oX* 
the  disturbance 

(Itf) 

.^-1 
*1  ^ 

(15) 

B 

Usve~veloclty  of 
the  disturbance 

(16)  a* 

O 

u, 

X 

Specific  heat  at 
constant  volume 

(17)  C 

V 

1 

n 

''v 

Specific  heat  at 
constant  pressure 

fifl  ^ r. 

-p 

y 

V 

°v 

Gas  constant  per  gram 

(19)  R"" 

K -1 

Acceleration  due  to 
gravity 

(20)  g 

1 

7 T '5 

V -X 

i. 

r 

o 

Mcxaenf  iT!  thickiuisj  of 
boundci'i'  layer 


UU 

9 - y p'^i  l^v)dx^ 

o 


(21) 
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( 


I 'I 


Dimensional  llon-dlmensional  Reference 

ouanti tlea  quantities  quahtlt ies 


Reync^lda  number 

( id ) Based  on  boundary 
layer  thickness  5 


i'dj)  Based  on  mcmentum 
thickness  0* 


R = 


R^ 

o 


7T 

‘W~ 


(24)  Based  on  distance  R 
from  leading  edge 


Mach  nvimber 

(25)  «,  = — ^ 


Proucle  number 

(26) 


Prandtl  number 


(27) 


k* 


Remarks 

(c  ) A bar  ( ) denotes  mean  value,  a dash  ( )*  denotes  fluctua- 

tiono  The  subscript  ( )^  refers  to  the  value  at  the  edge  of  the 

boundary  layer  at  the  particular  x-location  considered,  and  the 
subscript  ( refers  to  wall  value.  The  subscripts  ( )^,  and 

{ )j^  denote  the  reel  and  Imaginary  pa-^  ts  of  a quantity,  respectively. 


! 

I 
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(b)  me  non- dimensional  mean  quantities  UgCx,  .Xg 

^v3c*y)»  p(x,/),»»»ar©  repreaentad  as  flmeticne of  j alone,  that 
es  w(y),  v(y),  p(y ),p(y),- • •respectively.  This  corresponds  to  an 
approximation  depending  on  the  boundary-layer  neture  of  the  moan 
flow  (see  Section  l^.). 

(c)  The  quantities  C_,C  , yU,\,k,or,  y are  all  considered  to  be 
functions  of  T ^nly.  The  fluctuation  of  sny  such  quantity  Q 
Is  given  by  Q*  = ui*/ dif  T*. 

(d)  For  moderate  Mach  liurabers  and  rates  of  heat  transfer  C ,C  , y,cr 

Jr  • 

vai^  ail'y  sxlf^ly  in  the  boundary  layer.  They  are  therefore  re- 
garded as  constants  equal  to  their  free  stream  values  lii  the 
numerical  oaloulstions- 

(e)  The  reference  length  S is  the  boundary-layer  thickness 
at-:,  tbile  pcirtlcul&r  value  of  x considered.  Its  > definition  is 
somewhat  arbitrary,  and  may  be  taken  in  any  convenient  mannsr. 

In  the  present  numerical  computations  it  is  taken  to  be  the 
value  of  Xg  *t  which  * 0-999  (s®e  Appendix). 
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Table  2(b) 
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The  function  C = + B/A  for  zero  heat  transfer- 
er = 1,  « 1. 

0 0.5  0.9  1.3  1.7  2.i|-72 

C 0 O.I4.55  0.?'?9  0.9^^.1  1.006  0.885 

Table  5(b) 

The  function  C = V^7[T~c)*  S/A  for  various  temperature 
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(c)  The  I'unctton  w = l/2f * ( ^ ):  Jiee  prliii  "or  definition  of  f(  ^ ) 
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Purctlona  for  c8?.culatlon  of  curve  of  neu'sral  stability 
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